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ABSTRACT 


A  study  to  determine  the  effect  of  1)  data  power  spectrum 
and  2)  system  design  parameters  on  aliasing  and  data  interpolation 
error  has  been  conducted.  Results  are  applicable  to  time-division 
multiplexed  telemetry  systems.  In  particular,  the  results  apply  to 
the  commonly  used  PAM-FM,  PDM-FM  and  PCM-FM  telemetry  systems.  Data 
was  obtained  on  the  basis  of  combined  analysis  and  experiment  with 
emphasis  on  the  derivation  of  experimentally  proven  design  parameters. 
The  program  consisted  of  three  phases:  1)  A  survey  was  conducted  to 
determine  typical  telemetry  data  spectra  2)  Analysis  was  performed 
to  predict  accuracy  performance  of  practical  data  and  interpolation 
filters  using  selected  examples.  Included  is  a  more  general  analysis 
to  determine  optimum  filter  transfer  functions,  both  realizable  and 
non-reallzable,  3)  Experimental  tests  using  data  spectral  models 
were  performed  to  derive  proven  interpolation  accuracy  performance 
on  the  basis  of  practical  filters  and  design  parameters.  The  experi¬ 
mental  tests  utilized  a  PAM-FM  system  to  derive  an  accurate  measurement 
of  attainable  performance  and  to  indicate  when  other  equipment  factors 
begin  to  Impose  a  limit  on  data  recovery  accuracy. 
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PURPOSE 


OBJECTIVE 

The  objective  is  to  provide  telemetry  users  with  a  quanti¬ 
tative  set  of  experimentally  proven  operating  parameters  which  can  be 
utilized  to  avoid  a  significant  source  of  telemetry  data  error  (aliasing 
error)  in  time-division  telemetry  systems. 

TASKS 


Specific  efforts  undertaken  in  the  conduct  of  this  program  were: 

1 .  Determination  of  Typical  Data  Spectra 

A  study  was  conducted  of  actual  data  signals  to  determine 
typical  spectral  characteristics. 

2.  Extension  of  Past  Analyses 


An  extension  of  past  analyses  was  performed  to  better  pre¬ 
dict  aliasing  errors  in  actual  time-division  telemetry  systems. 

3.  Laboratory  Measurements 

Laboratory  measurements  were  conducted  on  actual  systems  to 
verify  and  extend  theoretical  unoer standing  of  the  error  mechanisms  in 
time-division  telemetering. 

Consideration  was  given  to  a  number  of  practical  design  para¬ 
meters  in  addition  to  the  specific  areas  of  effort  described  above.  Ex¬ 
amples  of  parameters  considered  include:  1)  data  sampling  rate,  2)  sampling 
aperture  duty  cycles,  3)  data  spectral  cutoff  rate,  4)  interpolation  fil¬ 
ter  cutoff  rate,  and  5)  the  interpolation  filter  bandwidth  normalized  with 
respect  to  the  data  bandwidth.  Of  particular  interest  are  the  test  results 
from  a  100%  duty  cycle  PAM  sampled  data  system.  This  system  employed 
essentially  instantaneous  data  samples  which  were  converted  to  100%  duty 
cycle  pulses.  The  original  data  was  reconstructed  from  a  multi-channel 
system  by  use  of  a  60%  aperture  sampler  followed  by  an  interpolation 
filter. 
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SECTION  1 


SUMMARY 


1.1  PROGRAM  OBJECTIVES 

A  study  program  was  conducted  to  1)  evaluate  the  effect  of 
errors  arising  from  "aliasing"  of  data  signal  spectra  by  a  sampling 
process  2)  investigate  the  relationship  between  the  aliasing  error  and 
distortion  error  produced  by  an  interpolation  filter^  and  3)  derive  a 
practical  understanding  of  and  a  proven  basis  for  predicting  the  error 
mechanisms  involved  in  the  acquisition,  transmission  and  reconstruction 
of  data  signals  by  means  of  a  sampled-data  telemetry  system. 

In  the  course  of  a  telemetry  system  study  recently  concluded 
(see  Reference  1)  it  was  found  that  little  information  of  a  practical 
design  nature  existed  which  took  into  account  the  aliasing  and  inter¬ 
polation  distortion  error  on  overall  system  accuracy.  The  complexity 
of  present  day  telemetry  instrumentation  with  the  inevitable  requirement 
for  handling  tens  or  hundreds  of  channels  of  information  makes  the  use 
of  time-division  multiplexing  almost  universal,  sometimes  in  conjunction 
with  frequency  multiplexing. 

When  the  proper  operating  parameters  are  not  certain  it  is 
common  practice  to  assign  sampling  rates  and  other  design  parameters 
thought  to  be  conservative  from  the  standpoint  of  minimizing  data  re¬ 
covery  errors.  However,  conservatism  carries  with  it  performance  penal¬ 
ties  involving  excessive  RF  spectrum  bandwidth  and  power  requirements. 

The  net  result  of  inadequate  knowledge  regarding  proper  selection  of  equip¬ 
ment  design  parameters  is  to  impose  inequitable  assignment  of  error  toler¬ 
ances  throughout  instrumentation  systems.  This  results  in  excessively 
stringent  specifications  on  equipment  in  some  areas  and  lack  of  adequate 
per formanci.'  in  others.  This,  in  turn,  leads  to  .sy.stem  designs  which  are 
excessively  sensitive  to  small  changes  in  equipment  performance. 
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Telemetering  systems  are  currently  designed  by  assuming  the 
data  source  has  a  known  bandwidth.  The  spectral  shape  of  the  data 
source  Is  rarely  taken  Into  account.  This  has  led  to  excessively  con¬ 
servative  designs  In  addition  to  the  simple  fact  that  accuracy  cannot 
be  adequately  predicted  even  on  the  basis  of  knowing  that  conservative 
design  practice  has  been  followed. 

Detail  consideration  of  the  source  of  design  difficulty  dis¬ 
cussed  above  led  to  the  program  described  by  phase  in  the  following 
paragraphs . 


1.2  DATA  SPECTRA  SURVEY 

Considerable  qualitative  information  was  found  as  discussed 
in  Section  2.  These  data  were  classified  into  five  spectral  types 
which  served  as  models  for  the  analytical  and  experimental  work.  A 
compilation  was  made  of  available  data  on  telemetry  data  spectra  from 
information  gained  by  means  of  a  literature  search,  telephone  inquiries, 
letter  inquiries,  and  personal  contacts. 


1.3  ANALYSIS 

Two  new  analytical  tasks  were  carried  out.  One  consisted  of 
a  general  PAM  system  analysis  including  such  practical  effects  as  finite 
sampling  width,  crosstalk,  and  noise.  It  specifies,  in  general  form,  the 
optimum  sampling  rate,  and  the  optimum  transfer  functions,  both  realiz¬ 
able  and  non-realizable ,  for  the  data  and  interpolation  filters.  The 
other  analytical  work  was  an  extension  of  past  analyses  to  include  the 
effects  of  practical  filter  transfer  functions,  and  the  effects  of 
interpolating  100%  duty  cycle  pulsetrains  as  well  as  impulse  or  short 
duty  cycle  pxilsetrains .  This  work  was  programmed  and  a  limited  amount 
of  computation  performed  with  an  IBM-709  computer.  The  analytical  work 
is  described  in  Section  3. 


1.4  EXPERIMENTAL  EVALUATION  PROCTIAM 

Laboratory  tests  were  performed  on  two  systems  which  give  in¬ 
formation  useful  in  design  of  a  time-division  multiplexed  telemeter  to 
handle  the  expected  data  types.  The  majority  of  tests  were  carried  out 
with  a  single  channel  P/\M  system  using  no  currier  since  this  simplified 
model  includes  the  parameters  of  importance  to  interpolation  error  which 
can  be  considered  to  be  composed  of  aliasing  and  distortion  errors.  At 
the  same  time  the  simplified  model  removes  error  sources  which  would 
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otherwise  obscure  the  relationships  of  Interest. 

Selected  tests  were  also  performed  on  an  8  channel  FAM-FM 
system  with  noise  added  so  that  all  basic  system  error  mechanisms  were 
present.  The  ..overall  system  error  was  determined  as  a  function  of  the 
signal-to-noise  ratio.  It  was  verified  that  overall  performance  of  such 
a  system  can  be  accurately  estimated  by  combining  the  interpolation 
error  from  single  channel  tests  with  the  pulse  sampling  error  resulting 
from  all  the  other  error  sources  in  the  system.  Pulse  sampling  error 
may  be  determined  for  the  three  common  time-division  systems  from  in¬ 
formation  obtained  in  Reference  1  or  a  similar  source  of  data.  The 
experimental  program  is  described  in  Section  4.  The  experimental  test 
setup  and  equipment  used  are  described  in  Section  5. 


SECTION  2 


DATA  SPECTRA  SURVEY 


A  survey  was  undertaken  to  determine  typical  spectral  character¬ 
istics  of  actual  data  signals  originating  from  information  sources. 
Knowledge  of  the  power  spectral  density  of  the  information  channels  is 
important  since  the  interpolation  error  is  strongly  dependent  on  the 
spectra  of  the  data  being  sampled.  A  literature  search,  inquiries  by 
telephone  and  letter,  and  personal  contacts  were  made  to  gather  the  data 
which  were  then  categorized  into  five  basic  spectral  types.  A  broad 
range  of  general  data  was  uncovered^  however,  specific  details  of  data 
characteristics  were  usually  missing.  This  led  to  recovery  of  a  large 
amount  of  information  not  amenable  to  reduction  to  specific  characteristics 
in  the  form  of  tabulated  results.  Therefore,  a  summary  of  the  most  useful 
references  and  sources  is  included  in  Appendix  3. 


2.1  BASIS  FOR  DATA  SPECTRA  SURVEY 

A  data  channel,  as  an  information  source,  is  completely  speci¬ 
fied  when  the  a  priori  symbol  probability  distribution  and  the  average 
rate  of  symbol  generation  is  defined.  Interpretation  of  the  term 
"symbo] "  is  determined  by  the  physical  representation  of  symbols  by  the 
channer  in  question.  These  symbols  may  be  discrete,  continuously  variable, 
or  a  composite  of  discrete  and  continuous  states  of  a  physical  parameter. 
The  corresponding  probability  distributions  associated  with  the  possible 
states  of  an  information  channel  may  be  discrete,  continuous  or  correspond 
to  a  mixed  process.  Selection  of  the  method  of  representing  a  symbol,  and 
the  mechanism  of  changing  from  one  symbol  to  another  coupled  with  the 
probability  distribution  of  the  symbols  permits  calculation  of  the  power 
spectral  density  of  the  information  source.  Since  the  power  spectral 
density  of  an  arbitrary  waveform  can  be  calculated  from  the  first  and 
second  moments  of  the  "symbol"  probability  distribution  it  is  evident 
th^  tjje  power  spectral  densi  tv  of  an  information  channel  is  specified 
hv  n,  .1^,  where  a,  are  the  first  and  second  moments  oi  the 

physical  parameter  representing  the  state  of  the  information  source, 

S(t,a)  the  function  of  time  acting  as  a  vehicle  for  carrying  the  symbol 
states  through  the  communication  link. 
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In  practice,  it  is  a  matter  of  definition  of  what  constitutes 
the  point  in  a  telemetry  link  where  information  is  generated.  For 
practical  instrumentation  this  corresponds  to  the  output  of  a  transducer 
of  a  specified  type.  Hence  the  function,  S(t,a)  is  well-defined  in  each 
individual  application. 

Conduct  of  the  data  survey  was  based  on  the  fact  that  transducer 
characteristics  are  strongly  established  and  therefore  should  generate 
spectra  typical  of  the  instrumentation  application. 

The  effect  on  power  spectral  density  by  the  information  source, 
is  reflected  by  the  distribution  of  energy  over  the  passband  of  the 
transducer.  In  those  cases  where  the  transducer  and  telemetry  link  do 
not  constitute  important  limits  on  spectral  distribution  of  the  observed 
source  it  is  reasonable  to  use  observed  instrumentation  data  to  determine 
the  spectra  generated  by  the  data  source. 


2.2  SURVEY  DATA  COLLECTION 

It  was  originally  thought  that  much  of  the  information  on  data 
spectra  could  be  obtained  through  ASTIA  in  report  form.  The  plan  was  to 
classify  the  data  source  as  to  transducer,  missile  type  and  transducer 
application,  thereby  obtaining  a  small  number  of  representative  spectra 
which  could  be  simulated  in  the  laboratory  as  typical  data  spectra.  This 
would  lead  to  realistic  system  operating  parameters  for  reconstructing 
the  transducer  outputs  to  within  a  specified  accuracy.  This,  in  turn, 
permits  recovery  of  the  data  symbol  "a"  to  within  a  specified  accuracy. 

Early  in  the  program  the  A.STIA  bibliographies  listed  in 
Appendix  4  were  ordered,  obtained  and  examined.  From  several  thousand 
bibliography  cards  examined  about  20  reports  were  estimated  to  have  the 
type  of  information  on  data  spectra  desired.  The  classified  reports  were 
not  ordered  because  of  security  restrictions  on  the  present  contract.  It 
is  interesting  to  note  that  the  majority  of  useful  reports  which  have  been 
examined  were  not  listed  in  the  ASTIA  bibliographies.  Several  reasons 
for  lack  of  information  available  through  ASTIA  are:  1)  it  is  difficult 
to  determine  the  proper  ASTIA  descriptors;  and  2)  the  reports  which  appear 
in  ASTIA  files  are  intended  to  be  of  general  usefulness  and  often  do  not 
include  detail  information  on  data  spectra. 
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Additional  sources  of  information  concerning  data  spectra  were 
needed  and  found.  Important  data  was  derived  from  instrumentation  estab¬ 
lished  by  Aeronutronic  on  the  Shillelagh  and  Blue  Scout  missile  programs. 
It  was  surmised  that  much  of  the  needed  information  existed  in  the  form 
of  internal  memoranda  and  reports  at  various  companies  and  military  test 
ranges.  A  survey  letter  was  prepared  and  mailed  directly  to  personnel 
in  22  companies  and  4  military  test  ranges  requesting  information  on 
data  spectra.  A  very  limited  number  of  replies  have  been  obtained  to 
dace. 


Several  problems  arise  in  evaluating  the  data  spectra  available 
in  the  reports  obtained.  These  are:  1)  determining  the  extent  to  which 
the  data  spectra  has  been  modified  by  the  transmission  i  item  used  to 
telemeter  the  data  when  little  or  no  mention  is  made  of  the  transmission 
system  used,  2)  establishing  a  means  of  comparison  of  data  spectra 
measured  by  different  techniques  which  are  not  usually  discussed  and 
3)  most  of  the  data  spectra  available  is  for  vibration  data  and  acoustic 
noise.  Information  about  other  flight  test  data  spectra  such  as: 
Acceleration,  Temperature,  Pressure,  Stress  and  Strain,  Velocity, 
Scientific  Measurements,  etc.  is  not  generally  available.  The  same, 
problems  were  encountered  in  another  study  discovered  during  this  survey 
carried  out  by  Fairchild  Engine  and  Aircraft  Corporation.  Their  work 
was  summarized  in  a  report  entitled  "Study  and  Evaluation  of  In-Flight 
Guided  Missile  Environmental  Data",  Final  Report,  15  December  1958, 
Fairchild  Engine  and  Aircraft  Corporation,  ASTIA  No,  AD  306  231. 

In  spite  of  the  limitations  in  much  of  the  data  collected,  it 
is  believed  this  is  a  representative  sample  of  that  available  and  it  is 
sufficiently  varied  to  permit  classification  into  typical  spectral  types 
permitting  the  simulation  of  data  spectra  for  the  laboratory  tests. 

One  problem  noted  above,  that  of  modification  of  the  data 
spectrum  by  the  transmission  system,  can  be  largely  removed  by  wire 
line  transmission  in  tie-down  tests.  One  test  report  referenced  in 
Appendix  3  describes  a  tie-down  test  with  measurements  of  vibration  and 
electrical  circuit  noise.  The  data  shown  indicates  a  significant  per¬ 
centage  of  energy  exists  up  to  and  beyond  6,000  cps  for  both  data  types. 
Therefore,  in  many  tests  it  must  be  assumed  that  the  high  frequency 
content  of  data  signals  such  as  vibration,  and  electrical  and  acoustical 
noise  are  limited  by  the  telemetry  system. 


Survey  results  have  been  reduced  in  usefulness  from  lack  of 
information  regarding  the  spectral  distribution  of  data  channels  requiring 
moderate  to  low  frequency  response.  This  is  perhaps  natural  but  prevents 
deducing  the  correctness  of  spectral  models  derived  from  consideration  of 
the  nature  of  the  data  sources  and  the  good  survey  results  involving 
mechanical  vibration,  acoustic  vibration  and  other  noise -like  wideband 
data.  In  general,  a  good  survey  sampling  was  obtained  for  wideband  data 
signals  which  currently  account  for  most  of  the  system  bandwidth  in 
telemetering  systems. 

Future  extension  of  knowledge  concerning  the  underlying  mech¬ 
anisms  generating  the  low  frequency  data  channels  would  permit  better 
telemetering  system  performance  and  at  the  same  time  gain  the  benefit 
of  significant  RF  bandwidth  reduction. 


2.3  DATA  SPECTRA 

It  was  found  that  regardless  of  source,  spectral  characteristics 
of  data  reported  or  predicted  can  be  categorized  into  five  types  of  spectra 
defined  by  a  set  of  parameters  with  variation  of  individual  parameters 
permitted  within  each  type.  The  classifications  do  not  exhaust  all 
possibilities;  however,  few  spectra  will  be  encountered  which  cannot  be 
realistically  represented  by  the  types  defined.  Thought  was  given  to 
categorizing  the  data  type  according  to  the  data  source,  e.g.  structures 
data  (vibration,  stress,  strain),  aerodynamic  (velocity,  acceleration, 
position),  environment,  subsystem  performance,  and  so  on.  However,  this 
method  of  classifying  did  not  appear  useful  to  the  study  or  control  of 
interpolation  errors  in  system  design  and  applications.  The  spectral 
shape  of  the  data  is  of  prime  importance  to  the  interpolation  error  and 
a  given  data  source  may  not  always  exhibit  the  same  spectrum.  For 
instance,  it  was  found  that  vibration  and  some  acoustic  noise  data 
spectra  can  be  represented  by  band-limited  white  noise  with  various 
cutoff  ratios  (db/oct.)  for  most  liquid  fuel  rockets  or  missiles,  while 
solid  and  some  liquid  fuel  rockets  have  vibration  and  acoustic  noise 
spectra  exhibiting  a  comb  spectral  structure. 

Practical  reasons  limit  consideration  to  and  definition  of  five 
general  spectral  types  discussed  below.  These  types  were  simulated  in  the 
laboratory  for  conducting  the  experimental  test  program. 
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The  spectral  types  are  considered  to  be  models  of  the  spectra 
generated  by  the  various  data  sources.  Only  positive  frequencies  are 
shown.  Classification  of  a  data  source  by  the  telemetry  design  and/or 
instrumentation  engineer  permits  entry  and  use  of  the  parametric  error 
curves  presented  in  Section  4. 


TYPE  I  -  Continuous  Spectrum 


f 


This  is  a  continuous  spectrum  flat  over  the  passband  defined. 

The  cutoff  characteristic  is  parameterized  and  permitted  to  attenuate 
at  6J  db/octave  where  J  is  a  positive  integer. 

Actual  data  sources  usually  have  a  low  frequency  cutoff.  This 
has  little  effect  on  the  accuracy  of  interpolation  or  reconstruction  of 
the  source  signal.  The  skirt  characteristic  beyond  the  cutoff  frequency 
largely  controls  the  seriousness  of  interpolation.  This  is  considered  the 
most  important  data  type  for  general  consideration  since  it  is  a  common 
one  encompassing  the  wideband  noise-like  data  sources.  This  type  repre¬ 
sents  much  of  the  vibration  and  acoustic  noise  data  observed  in  liquid 
fuel  missiles,  and  gust  turbulence  which  is  independent  of  the  vehicle 
considered.  The  skirt  cutoff  rate  should  be  estimated  by  the  instrumenta¬ 
tion  engineer  from  information  supplied  to  him  by  the  applicable  technical 
groups,  since  data  filtering  is  invariably  used  to  attenuate  transducer 
resonances  and  to  prevent  overmodulation.  The  chosen  parameter  J  then 
defines  the  cutoff  rate  of  the  source  spectrum  determined  by  the  applica¬ 
tion,  Thus,  two  parameters  ( )  define  the  spectral  shape.  Commonly 
encountered  ranges  for  these  parameters  are  o  -f^  ^  lo  ooo  •  |  -  J  £  $> 
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Type  II  -  Discrete  Spectra 


This  spectrum  is  made  up  of  a  sum  of  discrete  spectral  lines. 
Data  with  narrow  passband  (comb)  characteristics  may  be  represented 
by  such  a  discrete  spectrum.  Data  sources  of  this  type  typically  have 
less  than  12  spectral  components  of  importance  with  1  to  3  components 
usually  being  very  dominant.  This  type  is  typical  of  vibration  and  noise 
data  where  strong  resonances  are  excited.  In  the  majority  of  cases 
some  energy  is  observed  between  discrete  spectral  lines  where  structural 
or  acoustical  cavity  resonances  have  been  excited.  This  results  in 
addition  of  spectral  Type  I  with  Type  II  which  produces  a  mixture  of  a 
continuous  and  a  comb  spectrum.  This  results  in  a  Type  III  spectrum. 

It  was  predicted  and  confirmed  experimentally  that  wlien  the  energy  in 
the  continuous  spectrum  of  data  Type  III  drops  well  below  the  energy 
in  the  comb  spectrum  and  the  continuous  spectrum  cutoff  frequency  is 
not  above  the  highest  spectral  line  frequency  the  interpolation  error 
and  the  system  parameters  are  essentially  determined  by  the  comb  portion 
of  the  spectrum.  In  this  case  data  Type  III  can  be  simplified  to  Type  II. 

Type  III  -  Continuous  Plus  Discrete  Spectrum 


A  Type  III  spectrum  is  the  sum  of  a  continous  and  a  discrete 
spectrum.  The  spectral  model  is  defined  by  the  frequencies  of  the 
discrete  spectral  lines,  the  relative  power  in  each  line,  the  cutoff 
IroqLency  ■’!  the  continuous  spectrum  and  the  cutoff  characteristic. 

In  addition,  one  more  parameter  is  required,  the  ratio  of  power  in 
the  discrete  to  the  power  in  the  continuous  spectrum.  This  defines  the 
relative  levels  of  the  energy.  Restricting  attention  to  the  most 


9 


important  cases  where  the  dominant  spectral  lines  are  3  or  less  in 
number  permits  use  of  five  or  less  parameters  as  (td  i  ) 

to  specify  the  model  where  defined  as: 


fc  -  power  in  continuous  spectrum 
P;  =  power  in  i^^  spectral  line 

I 

The  set  of  five  parameters  define  the  model  completely  except  the 
total  power  of  the  spectrum  which  is  expressed  by  the  equation: 

-  Pc  +  IP.  --  Pc('  W.  +^3  ) 

a:  • 

This  is  an  arbitrary  parameter  which  does  not  effect  the  internal 
structure  of  the  model. 

Type  III  spectra  often  occur  where  a  broadband  noise-like 
energy  source  excites  system  resonances.  It  has  been  found  that  some 
engines  and  boosters,  and  gust  turbulence  are  essentially  white  noise 
energy  sources.  When  these  sources  excite  structural  resonances 
Type  III  spectra  are  seen  in  vibration  and  acceleration  measurements. 

A  Type  IV  spectrum  has  a  continuous  power  spectral  density 
with  a  shape  determined  by  two  critical  frequencies  as  illustrated 
in  the  diagram  below. 

Type  IV  -  Continuous  -  With  Step 
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It  does  not  appear  worthwhile  to  consider  and  define  spectral 
types  having  more  than  one  transition  from  one  constant  level  of  spectral 
power  density  to  another.  This  spectral  type  is  defined  by  specifying 
the  two  critical  frequencies, the  transition  attenuation  slope  the 
final  rate  of  attenuation  J2,  and  the  ratio  of  power  residing  in  the  two 
frequency  bands (  o  -  f*  )  and  (fa  —  )•  Thus,  five  parameters 

define  the  structure  of  the  spectrum. 

Type  IV  describes  actual  data  sources  having  much  of  their 
energy  in  the  lower  frequency  region  where  it  is  essentially  flat.  The 
other  portion  of  the  data  spectrum  of  lower  energy  but  still  having 
significant  energy  relative  to  that  in  the  low  frequency  region  may 
extend  up  1  or  more  octaves  beyond  the  first  critical  frequency.  This 
is  thought  to  be  a  good  model  for  many  data  types  which  are  essentially 
low  frequency  functions  but  vdiich  have  significant  higher  frequency 
perturbations.  Examples  are  accelerations  with  vibration  effects 
superimposed,  and  fuel  and  exhaust  pressures  with  added  engine  and  pump 
pulsations,  or  remote  signal  strength  measurements  with  fading  and  multi - 
path  effects.  This  can  also  be  considered  a  useful  model  for  continuous 
spectra  with  high  frequency  de-emphasis.  Some  acoustic  data  looks  like 
this, peaking  near  the  low  end  and  rolling  off  above. 

Type  V  -  Bandpass  Spectrum 


Type  V  are  continuous  bandpass  spectra  illustrated  in  the 
accompanying  figure.  This  type  may  be  considered  a  close  relative  of 
Type  IV  since  transition  from  the  first  spectral  region  (0  -  f ,  )  to  the 
second  (fj^  -  f2)  is  characterized  by  a  rise  in  spectral  power  density. 

Four  parameters  define  the  spectrum.  These  are  the  two 
critical  frequencies  (fj^  and  f2)  the  rising  slope  up  to  and  the 
falling  slope  J2  starting  at  f2  and  continuing  beyond.  Symmetric 
bandpass  characteristics  may  usually  be  assumed  which  makes  j J  j  =  I  J2 
and  the  parameters  are  reduced  to  three.  ^ 
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In  this  data  type  the  energy  is  concentrated  in  a  bandpass 
region  of  sometimes  only  one  octave  or  even  less.  This  data  type  is 
not  too  frequently  instrumented  so  that  less  emphasis  was  given  it. 
However,  some  experimental  measurements  were  made  since  it  is  an 
easily  distinguished  data  type.  Data  of  this  nature  was  found  in  the 
measurement  of  electrical  and  microphonic  noise  in  electronic  circuits. 
Some  of  the  data  found  were  actually  the  outputs  of  bandpass  filters 
in  the  electronics.  The  others  were  wider  band  but  still  of  a  bandpass 
nature. 
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SECTION  3 


ANALYSIS 


The  prime  objective  of  the  present  study  has  been  to  provide 
a  quantitative  set  of  experimentally  proven  operating  parameters  which 
can  be  utilized  to  avoid  a  significant  source  of  telemetry  data  error 
due  to  imperfect  interpolation  in  sampled  data  systems. 

The  number  of  methods  for  implementing  each  major  element  of 
a  sampled  data  system  are  usually  restricted  by  practical  equipment  de¬ 
sign  considerations.  Moreover,  it  is  evident  from  a  review  of  the  pos¬ 
sible  combinations  of  techniques  capable  of  being  used  in  a  single 
system  that  a  few  specific  configurations  must  be  assumed  for  analysis. 
These  configurations  can  then  be  optimized  in  performance  experimentally 
as  well  as  by  analysis.  Experimental  departure  from  theoretical  results, 
upon  evaluation,  displays  clearly  1)  the  limitations  of  equipment  and 
2)  limitations  of  the  assumptions  on  which  the  analysis  is  based.  Con¬ 
sequently,  analysis  conducted  in  this  program  served  to  supplement  and 
guide  the  experimental  program  as  well  as  to  point  out  new  possibilities 
worth  further  investigation. 

A  general  PAM  system  analysis  was  first  performed  to  serve  as 
the  basis  for  an  optimum  system  design.  The  following  summarizes  the 
results  of  the  PAM  analysis  which  appears  in  Appendix  1: 

a)  A  simplified  system  model  which  includes  the  effects  of 
finite  bandwidth,  finite  sampling  width,  aliasing,  cross¬ 
talk,  and  noise  is  developed  for  an  N-channel  system. 

b)  A  general  expression  for  the  total  mean  square  error  due 
to  all  the  above  effects  is  obtained. 
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c)  The  error  minimization  is  carried  out  fur  the  special  case 
in  which  all  channels  are  identical  and  a  signal  power 
constraint  is  imposed. 

d)  The  optimum  transfer  functions,  both  realizable  and  non- 
realizable,  are  obtained  for  the  data  filter  and  the  inter¬ 
polation  filter.  In  general,  these  transfer  functions  are 
bandlimited  and  can  only  be  approximated  in  the  practical 
case. 

e)  The  optimum  sampling  rate  is  specified  in  terms  of  the 
filter  transfer  functions. 

f)  The  optimum  sampling  delay  is  specified  in  terms  of  the 
video  filters  and  the  design  of  these  filters  to  reduce 
crosstalk  is  mentioned. 

From  this  analysis,  it  was  found  that  the  many  parameters  and  error 
mechanisms  of  a  multichannel  radio  transmission  system  obscures  the  funda¬ 
mental  aspects  of  interpolation  which  are  of  primary  interest  in  this 
study.  For  this  reason  the  general  PAM  analysis  appears  in  Appendix  1 
for  reference  and  certain  of  the  results  are  applied  to  the  analysis  of  a 
simplified  model  which  follows  in  this  section. 

Since  the  errors  due  to  interpolation  are  independent  of  the 
video  and  carrier  modulation  methods,  the  simplified  model  analyzed  in 
this  section  considers  a  single  channel  system  and  omits  the  carrier  por¬ 
tion  of  the  system.  The  results  obtained  here  for  interpolation  errors 
can  be  applied  to  any  sampled  data  system  (e.g.  PAM,  PDM,  PCM)  and  com¬ 
bined  with  the  other  error  sources  for  a  particular  system  (which  can  be 
obtained  from  the  results  of  the  Telemetry  System  Study,  Reference  1)  to 
obtain  overall  system  performance. 

The  simplified  system  model  is  described  in  Section  3.1. 

Section  3,2  discusses  the  analysis  of  fractional  duty  cycle  sampling; 
while  Section  3.3  presents  the  analysts  for  100%  duty  cycle  sampling.  The 
computer  program  utilized  to  evaluate  some  of  the  integrals  obtained  in  the 
analysis  is  discussed  in  Section  3.4  with  the  computer  results. 

"Alias-Free"  sampling  is  discussed  in  Reference  16.  However, 
the  method  of  implementing  the  sampling  technique,  and  interpolation  for 
practical  instrumentation  is  not  knowri  at  this  time.  Further  investigation 
in  this  direction  appears  worthwhile. 
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3.1  SYSTEM  MODEL 


The  general  system  shown  In  Figure  A-1-1  can  be  simplified  as 
shown  in  Figure  3-1.  This  also  represents  the  experimental  single  channel 
sampled  data  system. 


Data,  d(t)  - 

Sample 

& 

Noise,  N(t)  — 

- • 

Hold 

g(t) 


Interpolation 

Filter 

_ lii) _ 


f(t) 


Figure  3-1  Simplified  Single  Channel  Sampled  Data  System 


The  input  data  is  considered  to  be  the  input  to  the  sample  and 
hold  unit.  The  effects  of  data  filtering  are  included  in  the  analysis,  but 
are  not  evaluated  separately.  The  model  also  allows  noise  to  be  Introduced 
with  the  data. 

The  sampling  unit  gates  out  a  time  sample  t]^,  fixed  relative  to 
the  sampling  period  Tj  of  the  Incoming  signal  for  fractional  duty  cycle 
sampling;  the  hold  unit  permits  up  to  100%  duty  cycle,  l.e.  a  short  time 
sample  is  obtained  using  the  sampler;  this  sample  is  held  by  the  hold  unit 
until  the  next  sample  is  made.  Thus  a  100%  duty  cycle  PAM  wavetrain  can 
be  obtained. 

The  interpolation  filter  provides  an  output  which  approximates 
the  data  input  delayed  a  time  tQ. 

The  data  being  sampled  is  assumed  to  be  a  stationary,  random 
amplitude  function  of  time,  and  uncorrelated  with  the  noise. 

The  data  characteristics  considered  most  important  are  1)  the 
first  order  probability  distribution  and  2)  the  data  spectrum.  No  con¬ 
straint  Is  applied  to  the  data  spectrum,  i.e. ,  it  may  have  any  shape 
desired. 

The  symbols  introduced  above  and  throughout  this  Section  are 
defined  in  a  glossary  of  symbols  at  the  end  of  the  section. 
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3.2  INTERPOLATION  OF  FRACTIONAL  DUTY  CYCLE  PULSETRAIN 


The  following  determines  the  normalized  rms  error  between  the 
interpolation  filter  output  and  the  data  input  with  added  noise  inter¬ 
ference.  The  rms  error  is  determined  and  the  analysis  is  performed  in 
the  frequency  domain  with  power  spectra.  Analysis  is  based  on  the  data 
sampling  process  being  performed  by  a  multiplicative  aperture  although 
many  other  types  of  apertures  may  be  considered. 

3.2.1  Power  Spectrum  of  Sampler  Output 

The  multiplicative  sampling  aperture  is  defined  to  be  equiva¬ 
lent  to  taking  the  product  of  [d(t)  +  N(t)^  •  s(t)  =  g(t) ,  where  s(t) 

is  the  sampling  function  generating  the  apertures  shovn  in  Figure  2-3. 


The  power  spectrum  of  s(t)  is  determined  in  Appendix  2  and 
is  given  by  the  expression: 


note  that  if  A  =  JL  and  tj^ — *"0,  then  s(t)  becomes  an  ideal  sampling 
tl 

function,  i.e.,  a  train  of  delta  functions.  The  power  spectrum  for  this 
case  becomes 


The  power  spectrum  of  the  sampler  output  with  h(t)  *  d(t)  +  N(t) 
as  the  input  and  s(t)  as  the  sampling  function  with  A  =  1  is  derived  in 
Appendix  2.  The  result  is: 


16  - 


The  effect  of  the  sampling  function  is  exaggerated  In  the  figure  for 
Illustration.  Note  that  two  effects  may  be  noted:  1)  replicas  of  the 
original  data  and  noise  spectra  are  repeated  at  all  multiples  of  fg 
2)  the  replicas  of  the  noise  spectra  retain  their  additivity  with  the 
data  spectra.  Data  aliasing  errors  occur  due  to  interference  with  the 
desired  data  spectrum  by  these  data  and  noise  sideband  spectra  extending 
downward  and  upward  in  frequency  from  the  data  and  noise  spectral  replicas 

centered  about  +  fg,  +  2  fg . which  lie  within  the  passband  of  the 

data  interpolation  filter  centered  about  f=0.  The  aliasing  spectrum  is 
given  by: 


^  [fd  >  f'tO  3.4 


The  effects  of  linear  data  filtering  can  be  considered  if  Pc[(f) 
is  replaced  with 


where 

Pd(f)  is  the  data  input  power  spectrum. 

IdC^)!  Is  the  square  of  the  magnitude  of  the  data  filter  trans 
fer  function. 

Derivations  of  optimum  data  filters  are  given  in  Appendix  1  and 
by  Spllker,  reference  18. 


3.2.2  Power  Spectrum  of  Interpolation  Filter  Output 

The  power  spectrum  of  the  interpolation  filter  input  can  be 

written  as 
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where 


where 


x> 


p;  ’  "i 


»k « •  *0 

The  power  spectram  of  the  interpolation  filter  output  is 


P3CO 


=  Pu 


rcf) 

KO 


3-7 


+  P. 


N  jlCOp 


TC+^ 


is  the  square  of  the  magnitude  of  the  interpolation 
filter  transfer  function. 


3.2.3  Normalized  RMS  Error 


The  normalized  rms  interpolation  error  is  defined  as 


where 

€CM 


3-8 


3-9 


with  f(t)  and  d(t)  subject  to  the  constraint 


f\i)  -  cl^Ci-M  3-10 

The  constraint  normalizes  the  error  simultaneously  with 
respect  to  the  data  input  and  interpolated  output  in  terms  of  a 
reference  power  level.  This  constraint  is  imposed  automatically  on 
an  experimental  basis  to  prevent  normalization  errors  and  to  minimize 
the  interpolation  error  reading.  is  defined  by: 

_  00 

€(.i>  =  JPjCfjdP  3-11 

-00 
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In  addition,  the  following  relations  are  assumed: 


dc-t)  -  -  6Ct)  =  o 


3-12 


The  mean  values  are  assumed  to  be  zero  to  simplify  the  analysis; 
however,  this  does  not  limit  the  generality  of  results  obtained,  since  the 
assumption  that  3^)=oand  Tifc)  -O  amounts  to  a  d.c.  shift  in  data  signal 
reference  level. 

The  power  spectrum  of  the  error  is  determined  by  means  of  the 
Fourier  transform  of  60^)  which  can  be  written  as: 

«  -ju.t 

e  dt 

9  -(30 

3-13 

where 

F.(f)  is  the  Fourier  transform  of  the  data  input  and 
d 

N(f)  is  the  Fourier  transform  of  a  member  of  the  ensemble  of 
possible  input  noise  functions. 

The  above  can  be  re-written  in  the  following  form: 

>(«  -  C-T,')Fjff)|klc+> -c 

1- KI«)  ^(DiwhC^') 

This  equation  may  be  rewritten  in  terms  of  the  absolute  value  and 
phase  angle  of  the  data,  noise  and  interpolation  filter.  These  are 
denoted  by  ^  Cf"")  respectively. 


--  I 
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The  resulting  expression  Is: 


m 


3-14 


+  k 


rco 


f 

e 


relation: 


The  error  power  spectral  density  Is  obtained  by  means  of  the 


p(;f)=l,„  ± 

T-^0  ‘ 


F^(+) 
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where  the  bar  denotes  an  ensemble  average. 
Thus, 


:  Ci.'f  ljC«  |ij>»(%cn-iiT4t)*(irj 

+  k‘|r(f)f  { (y  ^c:  [stf-  \)  i- Pm«- ^,)1  A'it Ph») 

yy\o 

where  K  Is  a  gain  constant  associated  with  the  Interpolation  filter  to 
satisfy  the  constraint  that 

The  cross  correlation  between  data,  noise  and  sampling  functions 
are  all  assumed  to  be  zero. 


The  optimum  linear  interpolation  filter  transfer  function 
minimizing  0^^  can  be  determined  by  the  calculus  of  variations. 
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The  result  is: 


^ts)|^<^+)+phc+^'''(^')  Birv*" 

hTo 


3-17 


where 


A  derivation  of  this  result  is  performed  in  reference  18  and 
is  not  repeated  here.  It  should  be  noted  that  the  noise  is  not  sampled 
in  the  referenced  material  whereas  Equation  (3-17)  assumes  sampling  of 
the  noise  by  the  sample  and  hold  unit.  When  the  noise  is  added  following 
the  sampling  process  it  is  accounted  for  by  deleting  the  summation  of 
noise  terms  in  Equation  (3-17)  and  replacing  the  sum  with  the  unsampled 
noise  power  spectrum. 

The  effects  of  non-optimum  interpolation  filters  are  of  interest 
since  the  interpolation  filter  is  difficult  and  frequently  not  possible 
to  synthesize  for  non-trivial  data  spectra.  Therefore,  the  optimum 
solution  for  will  not  be  used.  However,  5]^*  is  minimized 

when  d^Cf)  ?  271^^0  .  The  phase  characteristic  of  a  well  designed  interpolation 
filter  is  approximately  linear  over  the  passband  accounting  for  the  average 
system  time  delay. 


The  mean  square  error  is: 

oe 

C'  r  J 

^  -0>3 


subject  to  the  constraint  that 


where 
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Let 


|ico 


Then 


2, 


Ki 


+<2; 


ai 


Substitution  of  (3-21),  (3-22)  and  (3-23)  into  (3-16) 


P,(f) 


I(f-) 


cif 


oC 


-L  K-f)  df 

-oO 


-)E 


2 

O 


3-22 

3-23 

3-24 

3-25 

gives : 

3-26 

3-27 
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where 


2  , 

“  normalized  mean  square  noise  error  ■  ^  V  ^  / 


■  normalized  mean  square  aliasing  error  *  K  i/%) 


<5- 


normalized  mean  square  distortion  error 

oO 


-  --(tD 

«*  ^  -aO 

The  error  terms  can  be  written  as  follows: 


Iff) 


E«- 


cJf 
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3-30 


'5‘  <■  1 


Using  Equation  (3-24)  the  total  normalized  mean  square  error 
can  be  written  as: 


E*-z^ 


I  - 


jpj(4)|l(^)ldf  I 


3-31 


Specific  examples  using  the  above  equations  are  considered  in 
Section  3.4.  The  results  obtained  by  using  numerical  integration  are 

also  presented  there. 
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3.3  INTERPOLATION  OF  1007.  DUTY  CYCLE  PAM  PULSETRAIN 


3.3.1  System  Model  Simplification 


The  analysis  to  follow  is  similar  to  that  given  for  fractional 
duty  cycle.  Therefore,  Figure  3-1  applies  and  some  steps  can  be  elimin¬ 
ated. 

Power  Spectrum  of  Sample  &  Hold  Output 

The  sampler  gates  a  small  time  sample  which  in  the  experiments  is 
restricted  to  47.  or  less  of  the  sampling  period.  This  sample  is  held 
until  the  next  sample  is  taken.  Thus  the  output  waveform  of  the  .sample 
and  hold  closely  approximates  a  1007.  duty  cycle  FAM  signal. 

The  power  spectrum  of  a  100%  duty  cycle  PAM  waveform  is  known 
to  be  (reference  2). 


where 


Tg  ■  sampling  period  (seconds)  =  1/fg 
fg  =  sampling  frequency 


a  >  random  variable,  the  pulse  amplitude 
a?  mean  square  value  of  a 
^  Ok  >  square  of  the  mean  value  of  a 

-e 

•  0  otherwise 


Note  that  tunctionally  the  sample  and  hold  output  spectrum  is 
invariant  with  respect  to  the  data  spectrum;  the  important  parameters 
are  the  mean  square  and  mean  values  of  the  data  which  vary  with  the  data 
amplitude  probability  distribution  and  the  sampling  period. 
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The  mean  Is  assumed  zero  to  simplify  the  power  spectrum  expression 
and  the  mean  square  value  can  be  replaced  with  the  more  familiar  symbol, 
z  o?-3[*=  5*  where  (T^  is  the  rms  value  of  the  input  data 

and  is  the  rms  value  of  the  noise. 

Thus 


FIGURE  3-4  POWER  SPECTRUM  AT  SAMPLER  OUTPUT  FOR  1007.  DUTY  CYCLE  PAM 


Note  that  the  spectrum  no  longer  consists  of  spectral  replicas 
of  the  noise  and  data  spectra  centered  at  0,  +  fg»  +  2  fg  -  -  - 
Instead,  the  shape  of  the  power  spectral  density  is  continuous  and 
independent  of  the  input  noise  and  data  spectra.  The  power  in  the 
continuous  spectrum  is  distributed  between  the  data  and  noise  in  the 
ratio 

The  power  spectrum  of  the  interpolation  filter  output  is: 


3-34 


where 


Che  square  of  the  magnitude  of  the  interpolation 

filter  transfer  function. 
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3.3.2  Normalized  RMS  Error 


The  normalized  rms  error  is  defined  in  Equations  3-8  through 
3-12.  The  power  spectrum  of  the  error  is: 


f^(f)  =  linn  ± 
T-9«6  T 


3-35 


|tC*)|  |F..(f;|cc5(6,(*)-2jH.) 

where  K  is  a  gain  constant  associated  with  the  interpolation  filter  to  satisfy 
the  constraint  that  -- j  _ > 

Using  the  calculus  of  variations  results  in  an  expression  for  the 
optimum  interpolation  filter  as  before. 


The  rms  error  is  minimized  when  ^ 

^  ^  L 


3-36 


&C^) 


where  the  above  solution  for  jmjj  is  the  optimum  interpolation  filter. 

In  order  to  determine  the  mean  square  error.  Equations  3-18 
through  3-20  will  be  used  along  with  the  following: 


3-37 
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The  mean  square  error  is; 

^  ?€  (f)  ^  3 

=  2«3"-  2(iJ  |iw  I  j|jf 

The  mean  square  normalized  error  is: 


Specific  examples  using  the  results  above  are  considered  in  the 
following  section.  The  results  obtained  by  numerical  integration  are  also 
presented. 


3. A  EXAMPLES  AND  COMPUTER  RESULTS 
3.4.1  Fractional  Duty  Cycle 


When  aliasing  and  distortion  are  the  only  error  mechanisms 
present  the  total  normalized  mean  square  error  becomes: 


E 


2 
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by  the  use  of  equation  3-27.  Using  the  value  of  K  determined  in  3-25 
and  the  definitions  following  3-27  results  in: 


E 


z 

D 
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where 


>  J  PjWi  J-f 

-  oo 

oO 

-  oO 

The  following  example  is  considered: 


Example 


Data  Spectrum 

Butterworth 

Interpolation  Filter 
Butterworth 

IC+) 

Kf) 

^  1 

,  . . 1 . . 

Vi  f  (♦/,;)»“ 

The  following  parameters  are  defined: 


T  -  f/fj 

«  -  f. 
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2  2  2 

The  parameters  defining  E.  ,  E  and  E  then  assume  the 
following  form: 


aO 


The  evaluation  of  the  total  mean  square  error,  in  terms  of 
the  component  errors  E^2  and  Eq^  Is  performed  by  use  of  a  programmed 
IBM-709  Computer.  The  normalized  errors  in  RMS  percent  of  full  scale 
are:  100  E,  100  E^  and  100  respectively. 

Results  of  computer  calculations  presenting  the  theoretical 
aliasing  and  total  rms  error  vs.  samples  per  cycle  of  the  data  cutoff 
frequency  is  shown  in  Figures  3-5  through  3-8.  The  numerical  limit  of 
integration  was  T  ■  f/fj  *  10.  This  results  in  an  overall  computation 
error  of  less  than  57,  for  the  case  ofJ*«K“lwlthTV  ranging  from  -2  to  +  2; 
N  >  1,  and  R  ■  0.02.  This  error  is  considerably  reduced  as  J  and  K  are  in¬ 
creased  since  less  energy  is  present  beyond  T  ■  10  for  the  faster  cutoff 
rates.  The  parameter  R  *  lift  was  fixed  at  0.02  since  the  experimental 
program  test  results  show  that  the  fractional  duty  cycle  is  not  a  sensitive 
parameter. 

Figure  3-5  shows  the  rms  error  vs.  samples  per  cycle  for  a 
Butterworth  interpolation  filter  with  K  *  6  and  fi/fd  •  2.  The  parameter 
J  defining  the  data  spectral  cutoff  rate  is  variable.  The  solid  curves 
correspond  to  the  aliasing  component  of  error,  and  the  dashed  curves 
correspond  to  the  total  error.  The  distortion  error  becomes  constant  at 
low  samples  per  cycle  and  is  the  limiting  error  source.  The  reduction  In 
distortion  error  is  quite  pronounced,  being  approximately  a  factor  of  5, 
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0  4  8  12  16  20 

SAMPLES  PER  CYCLE,  £s/fd 


FIG.  3-5  BUTTERWORTH  SHAPED  DATA,  -  0.02,  f^/fj  -  2,  K  «=  6 
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RMS  ERROR 


FIG.  3-6  BUTTERWORTH  SHAPED  DATA,  -  0.02,  *  2,  K  «  8 
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RMS  ERROR 


FIG.  3-7  BUTTERWORTH  SHAPED  DATA,  48  DB/OCT. ,  *  0.02,  f^/f^  -  4. 
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FIG.  3-8  BUTTERWORTH  SHAPED  DATA,  48  DB/OCT. ,  =  0.02,  =  1 
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as  the  data  cutoff  rate  Is  Increased  in  increments  of  12  db/oct. 


t 


In  Figure  3-6  the  cutoff  rate  on  the  interpolation  filter  has 
been  increased  by  12  db/oct.  The  aliasing  and  total  error  curves  agree 
closely  with  curves  which  are  approximately  comparable  in  Figure  3-5. 

The  rms  error  is  slightly  greater  in  Figure  3-6  because  of  increased 
distortion  error. 

Figure  3-7  is  similar  to  Figure  3-5  for  J  *  8  except  that  in 
Figure  3-7  f ■  4.  The  aliasing  error  has  increased  and  the  distortion 
error  has  been  reduced  significantly.  The  increase  in  aliasing  error  is 
due  to  the  wide  interpolation  filter  passing  the  sideband  (aliasing)  spectra 
at  low  samples  per  cycle.  The  curves  approach  each  other  as  fs/fj  is 
Increased.  The  lower  distortion  error  is  accounted  for  by  observing  that 
the  interpolation  filter  eliminates  less  data  energy  since  the  interpolation 
filter  cutoff  is  four  times  the  data  cutoff. 

Figure  3-8  can  be  compared  with  Figures  3-5  for  J  =  8.  In 
Figure  3-8  f^/f^j  “  1  as  compared  to  f =  2  in  Figure  3-5.  In  Figure  3-8 
the  aliasing  error  has  been  reduced  and  the  distortion  error  has  increased. 
The  reduction  in  aliasing  error  is  attributed  to  the  fact  that  the  inter¬ 
polation  filter  rejects  more  of  the  energy  in  the  aliasing  spectra  for 
fj^/fj  =  1  than  for  f^/fj  =  2.  The  increased  distortion  is  due  to  the 
interpolation  filter  rejecting  more  of  the  data  energy.  These  calculations 
have  shown  good  agreement  with  the  experimental  results  reported  in 
Section  4. 


3.4.2  1007.  Duty  Cycle  Sampling 


By  use  of  Equation  3-39,  the  expression  for  total  mean  square 
error  simplifies  to;  oO. 


1 


1  if 

"TNtT  j 
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when  noise  interference  is  reduced  to  zero. 


Below  is  given  the  expressions  comparable  to  the  fractional  duty 
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cycle  example  presented  In  Paragraph  3.4.1  when  a  Butterworth  data  and 
Butterworth  interpolation  filter  spectral  characteristic  is  assumed. 


The  following  parameters  are  defined: 

Q  =  fs/£d  T  =  f/fd 

N  =  fi/fd 


o4 


These  parameters  are  used  in  programming  of  the  IBM-709  computer. 
The  normalized  rms  error  in  percent  is  100  E. 
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Several  calculations  were  made  using  an  upper  limit  of  T  «  10 
for  integration.  The  parameter  values  selected  were  N  «  2,  J  ■  4,  and 
K  ■  6.  The  calculated  results  differed  from  those  obtained  experimentally. 
The  discrepancy  has  not  been  completely  resolved  but  is  believed  to  be  du.e 
to  an  error  in  progranming.  Since  the  calculations  performed  for  the 
fractional  duty  cycle  case  agree  well  with  the  experimental  results,  the 
experimental  results  for  100^  duty  cycle  have  been  verified. 

GLOSSARY  OF  SYMBOLS 


Symbol 


det) 


H*) 


fd 


9C^) 


ICf) 

J 


K 

m) 


Definition 

Data  time  function 

Normalized  rms  error 

Normalized  mean  square  aliasing  error 

Normalized  mean  square  distortion  error 

Normalized  mean  square  noise  error 

Interpolation  filter  output  time  function 

Data  cutoff  frequency  (3  db  point) 

Interpolation  filter  cutoff  frequency  (3  db  point) 
Sampling  frequency 

Fourier  transform  of  (  )  time  function 
Sampler  output  time  function 
Interpolation  filter  transfer  function 
Data  cutoff  rate  (J  x  6  db/oct.) 

Interpolation  filter  cutoff  rate  (K  x  6  db/oct.) 
Noise  time  function 
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GLOSSARY  OF  SYMBOLS  (Continued) 


Symbol 


to 


t, 


Ta 

€C-t) 

(Te 


^6 

'^f 


Definition 

Power  spectrum  of  (  )  time  function 
Sampling  time  function 
System  time  delay 

Sample  pulsewidth  (fractional  duty  cycle) 

Sample  pulse  period 

Error  time  function 

RMS  error 

RMS  data 

RMS  interpolation  filter  output 
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SECTION  4 


EXPERIMENTAL  EVALUATION  PROGRAM 


4. 1  INTRODUCTION 

This  section  reports  the  results  of  the  experimental  evaluation 
program  to  provide  a  quantitative  set  of  experimentally  proven  interpola¬ 
tion  design  parameters.  These  parameters  can  be  utilized  to  avoid  a 
significant  source  of  data  recovery  error  in  existing  telemetry  systems. 
The  data  reported  was  obtained  by  setting  up  and  operating  the  equipment 
described  in  Section  5. 

Two  sampled  data  systems,  single  and  multi-channel,  were 
investigated.  This  permitted  verifying  and  supplementing  theoretical 
understanding  of  aliasing  errors  and  its  relationship  to  the  other 
components  of  error  in  time-division  multiplexed  telemetry.  These  two 
systems  were  set  up  and  tested  under  actual  operating  conditions  with 
simulated  data  signals  and  noise  interference  introduced. 

The  laboratory  results  for  the  two  systems  are  presented  in 
Sections  4.3  and  4.4  following  the  methods  and  procedures  used  in  the 
experimental  program  discussed  in  Section  4.2. 


4.2  METHODS  AND  PROCEDURES 

The  test  program  consisted  of  setting  up  the  single  and  multi¬ 
channel  sampled  data  systems,  supply  random  input  signals  typifying  actual 
data  spectra,  and  to  measure  the  system's  performance  with  and  without 
interference  introduced.  System  performance  was  measured  by  means  of  an 
accurate  error  comparator  which  repetitively  sampled  the  input  and  output 
of  the  system  at  corresponding  points  and  compared  the  samples  for  error. 

The  adjustments  of  the  error  comparator  enable  the  operator  to  remove  the 
effects  of  time  delay,  as  well  as  amplitude  and  bias  changes  within  the  system. 


The  system  design  parameters  of  prime  importance  are:  1)  the 
ratio  of  samples  per  second  to  the  data  bandwidth,  fg/f^,  samples  per 
cycle,  2)  interpolation  filter  and  data  cutoff  rates,  K  and  J  in 
units  of  6  db/oct.,  3)  the  ratio  of  interpolation  filter  bandwidth  to 
data  bandwidth,  and  4)  normalized  sample  pulsewidth,  tj^fj^.  These 

symbols  were  used  in  Section  3.  Particular  emphasis  was  placed  on 
Investigating  the  effects  of  these  parameters  using  practical  interpolation 
filters  under  actual  operating  conditions. 

The  single  channel  tests  were  performed  without  a  carrier  system 
or  added  noise  interference  to  permit  thorough  understanding  of  the  effects 
the  various  parameters  have  on  interpolation  errors. 

A  carrier  system  was  used  for  the  multi-channel  tests  with  added 
noise  interference.  The  "carrier"  system  utilized  an  FM-FM  40  KC  sub- 
carrier  oscillator  and  discriminator  equipment  followed  by  a  Gaussian 
video  filter.  This  is  the  same  equipment  utilized  for  PAM  measurements 
in  the  telemetry  study  program  (Reference  1).  It  was  adapted  to  the  PAM 
pulse  rate  used  (4  KC)  in  the  current  study. 

Nearly  all  of  the  data  described  in  the  following  sections  have 
been  normalized  to  dimensionless  ratios.  The  normalized  data  permit  the 
results  shown  to  be  scaled  to  predict  performance  for  a  system  of  different 
data  rate  and  bandwidth.  In  this  connection  it  should  be  observed  that 
results  may  be  scaled  without  regard  to  carrier  frequency.  Use  of  most 
of  the  normalized  parameters  are  obvious  except  perhaps  the  normalized 
signal-to-noise  parameter,  S/kjFg^,  used  for  the  multi-channel  system. 

The  relations  between  S/k]^Fg%  and  the  more  commonly  used  S/kj^  and  S/N 
are  given  by: 
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where 


S  =  rms  carrier  voltage 

k]_  “  rms  noise  voltage  per  unit  root  bandwidth 

Fg  ■  system  sampling  frequency,  equivalent  to  the  pulse  rate 

B  ■■  system  design  bandwidth 

N  ■■  rms  noise  voltage  in  bandwidth  B 


4.3  SINGLE  CHANNEL  SAMPLED  DATA  SYSTEM 

4.3.1  System  Description 

The  general  system  block  diagram  is  given  in  Section  3, 

Figure  5-1  along  with  a  general  system  and  detailed  equipment  description. 
Briefly  the  system  consists  of  a  data  source  followed  by  a  filter  denoted 
as  a  data  filter  to  provide  spectral  shaping  in  conformance  with  the 
desired  spectral  data  model.  The  data  filter  output  is  sampled  and  the 
samples  are  supplied  to  an  interpolation  filter.  The  rms  error  is  obtained 
by  comparing  the  interpolation  filter  output  with  the  spectral  shaping 
filter  output.  It  must  be  emphasized  that  for  most  tests  the  data  filter 
was  used  for  spectral  shaping  of  broadband  white  noise  with  a  specified 
amplitude  probability  distribution;  its  output  was  considered  as  the  data 
source.  This  is  equivalent  to  data  filtering  of  an  actual  data  signal 
only  when  the  signal  has  a  spectral  characteristic  which  is  flat  over 
the  passband  of  the  filter  and  several  octaves  beyond. 

The  system  parameters:  1)  samples  per  cycle  2)  sampling  duty 
cycle  3)  interpolation  filter  data  cutoff  rate  and  4)  the  ratio  of 
interpolation  filter  bandwidth  to  data  bandwidth  -  were  investigated  for 
selected  data  spectral  models  and  selected  interpolation  filter  types. 

4 . 3. 2  Test  Results 


At  least  one  simulation  of  each  of  the  five  data  types  described 
in  Section  2  was  made  and  performance  measured  for  various  sytem  parameters 
during  the  test  program.  The  test  results  are  discussed  in  the  order  in 
which  they  were  performed  as  listed  on  the  following  page. 


-  41 


Data  Spectrum 


Basic  Type 


1)  Single  slnewave  II 

2)  Butterworth  shaped  data  I 

a)  Various  cutoff  rates,  24,30,42,48 

b)  2-amplltude  probabilities,  Gaussian 
and  Equl -probable 

3)  One  Slnewave  plus  Butterworth  shaped  data  III 

4)  Three  Slnewaves  plus  Butterworth  shaped  data  II  &  III 

5)  Continuous  with  step  IV 

6)  Bandpass  data  V 


Spectral  measurements  and  waveform  photographs  taken  at 
various  points  In  the  system  are  presented  in  Section  4.3.3. 

Slnewave  Test  Results 


The  percentage  rms  error  versus  samples  per  cycle  of  a  fixed 
slnewave  frequency  Is  shown  in  Figure  4-1.  These  tests  served  also  to 
check  proper  system  operation  and  to  determine  the  dynamic  residual  error. 

The  curve  for  “  1  results  from  using  a  100  cps  slnewave 

and  a  data  sample  pulsewidth  into  the  interpolation  filter  of  20  /I  sec. 

A  100  cps,  36  db/octave  Butterworth  filter  was  used  for  interpolation. 

The  curve  for  f^/fj  “  4  results  from  using  a  50  cps  slnewave, 
a  data  sample  pulsewidth  of  20  yttsec.  converted  to  a  100%  duty  cycle 
pulsetrain  to  the  interpolation  filter.  The  filter  is  a  modified  200  cps 
Butterworth  filter  with  a  400  cps  trap  and  a  slope  of  60  db/octave 
around  300  cps.  This  filter  will  be  denoted  as  BWT. 

For  the  parameter  greater  than  8,  the  two  curves  are  nearly 

coincident  and  the  minimum  residual  system  error,  achieved  for  20  samples 
per  cycle,  is  the  same,  0.14%.  However,  for  samples  per  cycle  less  than  8, 
the  rms  error  for  f^/f^j  =  4  is  higher  than  that  for  fp/f^j  =  !•  The  major 
reason  for  this  difference  is  that  fg  *  for  the  poorer  curve.  Conse¬ 
quently,  the  interpolation  filter  passes  undesired  spectral  components. 
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FIG.  4-1  SINEWAVE  TESTS 
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This  indicates  that  for  practical  filters  fg  should  be  greater  than  fj^, 
and  later  results  show  that  fg  ^  ^  ^i  good  practice. 

Sampling  frequency  feed  through  was  observed  at  low  values  of 
fg/fd  the  fractional  duty  cycle  test.  Therefore,  the  BWT  filter  was 
designed  to  trap  the  sampling  frequency  at  frequencies  equal  to  or  greater 
than  400  cps.  The  phase  and  amplitude  characteristics  of  this  filter  are 
considerably  different  from  the  Butterworth  filter.  The  amplitude 
characteristic  approximates  a  rectangular  cutoff  more  closely  than  the 
Butterworth  filter.  However,  the  phase  characteristic  is  quite  non-linear 
close  to  the  corner  frequency.  The  amplitude  response  is  shown  in  Section 
5  for  all  filters  used. 

Most  of  the  following  tests  were  run  with  the  BWT  filter.  How¬ 
ever,  distortion  errors  were  found  to  be  the  limiting  factor  in  achieving 
low  rms  error  with  this  filter  due  to  its  poorer  phase  response.  There¬ 
fore,  Gaussian  and  Butterworth  filters  were  also  investigated  and  found 
to  give  lower  rms  errors  above  values  of  4  to  6  for  fg/fd  where  interpola¬ 
tion  filter  distortion  has  a  more  pronounced  effect  than  aliasing  on  the 
interpolation  errors. 

Butterworth  Shaped  Data  Test  Results 

The  analysis  indicates  that  the  aliasing  error  depends  only  on 
the  first  and  second  moments  of  the  amplitude  source  and  not  the  particu¬ 
lar  amplitude  probability  distribution  of  the  data.  Two  noise  sources  were 
used  to  verify  this  experimentally,  equi-probable  and  a  Gaussian  amplitude 
probability  distribution.  These  are  denoted  by  EPA  and  GAP  respectively. 
Identical  tests  were  conducted  with  both  noise  sources  being  spectrally 
shaped  by  a  Butterworth  filter.  The  results  for  these  tests  are  discussed 
below. 

Butterworth  Shaped  EPA  Data 

Figures  4-2  through  4-5  present  the  rms  error  versus  samples 
per  cycle  of  the  data  cutoff  frequency  for  data  sources  with  equi-probable 
amplitude  distributions  (EPA  Data  Sources) .  'Hie  parameters  on  each  curve 
are  1)  the  ratio  of  interpolation  filter  bandwidth  to  data  cutoff  fre¬ 
quency,  fi/fd>  2)  the  ratio  of  sampling  frequency  to  interpolation 
filter  bandwidth,  fg/f^.  The  parameter  wliich  varies  between  figures  is  the 
sample  pulsewidth. 

The  tests  were  performed  using  the  BWT  interpolation  filter, 
holding  the  data  cutoff  frequency  fixed,  and  varying  the  sampling  frequency. 
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RMS  ERROR 


FIG.  4-2  BUTTERWORTH  SHAPED  DATA,  EQUI-PROBABLE  AMPLITUDE,  =  0.01 
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RMS  ERROR 


FIG.  4-4  BUTTERWORTH  SHAPED  DATA,  EQUI-PROBABLE  AMPLITUDE,  -  0.1 
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FIG.  4-5  BUTTERWORTH  SHAPED  DATA,  EQUI -PROBABLE  AMPLITUDE,  100%  DUTY  CYCLE 
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It  is  apparent  from  these  curves  that  the  sample  pulsewidth  has 
relatively  little  effect  on  the  rms  error.  This  is  illustrated  in  Figure  4-6 
where  the  rms  error  versus  pulsewidth  is  shown  for  fixed  and  If 

the  pulsewidth  times  the  interpolation  filter  bandwidth,  tj^fj^,  is  less  than 
0.5,  the  interpolation  filter  output  is  quite  small  relative  to  the  data  input 
and  considerable  amplification  of  the  interpolation  filter  output  is  needed. 
This  amplification  Increases  the  system  residual  error  because  of  increased 
circuit  non-linearities  and  noise.  Amplifeation  of  the  interpolation  filter 
output  can  be  eliminated  if  100%  duty  cycle  is  used  since  this  waveform 
supplies  adequate  power  to  the  interpolation  filter. 

As  the  ratio  of  f^/f^  is  increased,  the  rms  error  curve  shifts  to 

the  right  and  down.  This  is  because  the  knee  of  each  curve  corresponds  to 

the  point  at  which  the  sampling  frequency  equals  the  tuned  frequency  of  the 
trap  in  the  BWT  interpolation  filter.  The  aliasing  spectrum  is  less 
attenuated  by  the  filter  increasing  the  rms  error  at  reduced  samples  per 
cycle.  For  fg/fj  samples  per  cycle  above  that  at  the  knee  the  aliasing 
spectrum  is  greatly  attenuated  and  the  curve  is  asymptotic  to  the  distortion 
error  plus  system  residual  errors.  The^ distortion  error  results  from  the 
non-linearities  in  the  amplitude  and  phase  characteristics  of  the  interpolation 
filter.  Lower  distortion  error  results  as  the  ratio  of  f^/fj  is  increased. 

Figures  4-7  and  4-8  permit  selection  of  the  minimum  sampling 

frequency  and  corresponding  interpolation  bandwidth  to  be  chosen  for  a 
specified  rms  interpolation  error  and  data  bandwidth.  Since  the  data  is 
insensitive  to  pulsewidth  only  one  curve  is  shown  for  fractional  duty  cycle. 

For  fixed  rms  error  the  minimum  point  on  the  curve  corresponds  to 
the  proper  sampling  frequency  and  corresponding  interpolation  filter  band¬ 
width  for  the  given  interpolation  filter  design.  To  the  left  of  the 
minimum  a  sharp  increase  in  samples  per  cycle  is  required  to  maintain  a 
constant  rms  error;  since  the  narrower  interpolation  filter  bandwidth 
greatly  increases  the  distortion  error.  To  the  right  of  the  minimum  the 
samples  per  cycle  must  increase  with  the  wider  interpolation  filter  band¬ 
width  to  reduce  aliasing  error. 

The  minimum  samples  per  cycle  is  slightly  less  for  100%  duty 
cycle  than  for  fractional  duty  cycle  Interpolation  for  the  same  rms  error. 

The  curves  shift  down  and  to  the  left  as  the  rms  increases, 
decreasing  the  samples  per  cycle  and  interpolation  filter  bandwidth. 
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FIG.  4-6  RMS  ERROR  VS.  NORMALIZED  PULSEWIDTH 
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no.  4-7  SAMPLES  PER  CYCLE  VS.  NORMALIZED  ZNTERPOLATIOM  FILTER  BANDWIDTH,  -  0.02 
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NORMALIZED  INTERPOLATION  FILTER  BANDWIDTH,  ft/^d 


FIG.  4-8  SAMPLES  PER  CYCLE  VS.  NORMALIZED  INTERPOLATION 
FILTER  BANDWIDTH,  lOOX  DOTY  CYCLE 
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NORMALIZED  INTERPOLATION  FILTER  BANDWIDTH,  fi/fd 


Butterworth  Shape  GAP  Data 


Figures  4<-9  through  4-12  present  the  rms  error  versus  samples  per 
cycle  for  data  sources  with  Butterworth  spectral  characteristics  and  a 
Gaussian  amplitude  probability  (GAP)  distribution.  The  same  BVTT  interpola¬ 
tion  filter  is  used  as  for  Figures  4-2  through  4-5.  A  data  cutoff  rate  of 
42  db/oct.  was  used  to  obtain  the  experimental  results  for  the  parametric 
curves  given  by  fi/f(j  ■  1  and  4.  The  data  cutoff  rate,  was  48  db/oct. 
for  the  fl/fd  ■  2  parametric  curves.  The  different  data  cutoff  rates 
resulted  from  the  cutoff  characteristics  of  the  noise  source  described  in 
Section  4.3.3. 

I  The  results  agree  closely  with  those  for  equi-probable  amplitude 
distribution  in  spite  of  the  significant  difference  between  these  two 
distributions.  An  apparent  exception  occurs  for  the  curves  with  fi/fj  “  2; 
however,  the  cutoff  rates  differ  and  therefore  are  not  comparable. 

Comparison  of  data  for  the  two  distributions  is  presented  in 
Figures  4-13  through  4-15.  These  figures  are  replots  of  Figures  4-2  through 
4-5  and  4-9  through  4-12  for  the  various  ratios  of  interpolation  bandwidth 
to  data  bandwidth.  This  combines  the  two  distributions  on  one  set  of  scales 
for  comparison.  These  replots  show  the  close  agreement  between  the  equal 
amplitude  probability  distribution  and  the  Gaussian  probability  distribution 
noise  for  fi/fj  “  1  and  4.  Therefore,  it  is  concluded  that  the  test  results 
are  not  significantly  affected  by  the  amplitude  distribution.  It  should  be 
noted  that  distortion  error  resulting  from  amplitude  non-linearity  in  the 
system  would  be  affected  by  the  amplitude  probability  distribution. 

For  fl/fd  “  2  the  test  results  using  Gaussian  data  achieved  a 
lower  minimum  rms  error  than  for  equal  amplitude  data  due  to  the  sharper 
cutoff  rate  of  data  signal.  For  the  Gaussian  data  source  as  for  equi- 
probable  amplitude  data,  the  sample  pulsewidth  is  not  a  sensitive 
parameter.  The  close  grouping  of  the  curves  for  constant  f^/fd  indicates 
this  insensitivity. 

Figure  4-16  serves  to  show:  1)  sample  pulsewidth  did  not 
significantly  affect  interpolation  accuracy  over  a  large  range  of 
parameter  values,  2)  the  probability  distribution  of  the  data  does  not 
affect  interpolation  accuracy. 

Figure  4-16  was  derived  from  Figures  4-13  through  4-15  by  averaging 
with 'respect  to  pulsewidth  for  constant  fi/fj.  The  averaging  included  both 
equi-probable  and  Gaussian  probability  distribution  data  with  f^/fd  ■  1  and  4. 
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RMS  ERROR 


FIG.  4-12  BUTTERWORTH  SHAPED  DATA,  GAUSSIAN  PROBABILITY  DISTRIBUTION,  1007.  DUTY  CYCLE 
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RMS  ERROR 


FIG.  4-13  BUTTERWORTH  SHAPED  DATA,  fi/fj-l,  J-7 
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RMS  ERROR  VS.  SAMPLES/CYCLE 
INTERPOLATION  FILTER:  BWT 

fi/fj-2 

DATA; 

PROBABILITY  DISTRIBUTION 

_ EQUI -PROBABLE,  J-7 

-  GAUSSIAN,  J“8 


tjf  -0.02 


tjfi=0.01 


V  V 

\  N 


t  fi=0.1 


SAMPLES  PER  CYCLE,  fg/f^ 


FIG.  A-14  BUTTERWORTH  SHAPED  DATA,  f^/f^  «  2 
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FIG.  4-15  rUTTERWORTi  SHAPED  DATA,  q/f^=4,  J=/ 
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RMS  ERROR 


FIG.  4-16  BUTTERWORTH  SHAPED  DATA  CURVES  AVERAGED  WITH 

respe';t  to  sampij:  pulsewidth 
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Separate  averaging  was  performed  for  f^/fs  ■  2,  since  the  data  source 
spectral  cutoff  rates  differ.  Here  the  higher  cutoff  rate  reaches  a 
lower  minimum  error  as  expected.  For  fi/f^  >  4  an  average  was  included 
in  the  figure  for  two  cases:  1)  fractional  and  2)  100%  duty  cycle 
because  a  significant  difference  was  present.  This  can  be  attributed  to 
the  difference  between  the  fractional  duty  cycle  sampling  and  the  100%  duty 
cycle  sampling  techniques. 

Figure  4-17  gives  the  rms  error  versus  samples  per  cycle  for  a 
data  cutoff  rate  of  30  db/oct.  (J“5)  with  the  BWT  filter,  and  the  parameters 
fj^/fj  =  2  and  tjf^  =  0.02.  The  results  for  the  100%  duty  cycle  pulsetraln 
are  plotted  in  the  same  figure  for  comparison.  When  these  curves  are  com¬ 
pared  with  Figures  4-10  and  4-12,  it  is  apparent  that  the  data  cutoff  rate 
Is  a  sensitive  parameter  affecting  both  aliasing  and  distortion  errors. 

In  Figure  4-17  the  rms  error  increased  by  a  factor  of  2.8  over  that  for 
Figure  4-12  for  100%  duty  cycle  sampling  for  the  parameter  fs/fj  greater 
than  10,  while  the  Increase  is  a  factor  of  1.6  for  tj^f^  =  0.02  over  that 
shown  in  Figure  4-10.  Distortion  error  takes  effect  at  small  values  of 
fg/fjj  in  Figure  4-17  than  for  Figure  4-10  or  4-12.  This  causes  the 
sharper  corner  at  about  5  samples  per  cycle. 

Figures  4-18  and  4-19  give  the  rms  error  versus  samples  per  cycle 
for  24  and  36  db/oct.  Butterworth  interpolation  filters.  The  two  sets  of 
curves  in  each  figure  are  for  data  cutoff  rates  of  30  and  48  db/octave. 
Interpolation  of  tjfi  “  0.02  and  100%  duty  cycle  pulsetrains  are  shown  for 
each  data  cutoff  rate.  The  100%  duty  cycle  pulsetrain  gives  the  best 
results  for  aliasing  with  the  24  db/oct.  interpolation  filter  while 
tj^fj^  ■  0.02  is  best  for  the  36  db/oct.  interpolation  filter.  The  difference 
in  observed  distortion  error  is  not  significant.  It  can  be  seen  that  the 
greater  data  cutoff  rate  allows  a  smaller  rms  error  for  a  fixed  samples  per 
cycle  as  expected.  Comparison  of  the  three  pairs  of  curves  for  30  db/oct. 
data  cutoff  shows  that  the  error  is  fairly  insensitive  to  the  interpolation 
filter  response.  The  BWT  filter  performs  best  at  low  samples  per  cycle, 
fs/fj,  due  to  its  greater  selectivity  as  noted  earlier.  The  36  db/oct. 
Butterworth  response  offers  better  performance  at  high  values  of  fg/fjj  due 
to  reduced  distortion  error.  Performance  is  still  quite  good  with  the 
24  db/oct.  filter  and  this  would  be  a  useful  design  for  low  frequency  inter¬ 
polation  when  an  RC  filter  is  desired. 

Figure  4-20  presents  the  rms  error  vs.  samples  per  cycle  for  a 
24  db/oct.  Gaussian  interpolation  filter  with  100%  duty  cycle  pulsetrain. 

The  data  cutoff  rate  is  48  db/oct.  and  f ^/fj  =  1  and  2.  The  Gaussian  filter 
provides  poorer  interpolation  than  any  of  the  other  three  filters  suggesting 
that  amplitude  response  is  of  more  importance  than  phase  linearity. 
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FIG.  4-17  BUTTERWORTH  SHAPED  DATA,  GAUSSIAN  PROBABILITY  DISTRIBUTION,  J=--5 
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FIG.  4-18  •  BUTTERWORTH  SHAPED  DATA,  GAUSSIAN  PROBABILITY  DISTRIBUTION, f^/fj.2  K-4 
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FIG.  4-19  BUTTERWORTH  SHAPED  DATA,  GAUSSIAN  AND  EQUI-PROBABLE  AMPLITUDE 
PROBABILITY  DISTRIBUTION,  f^/fj  =2  K  =  6 
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FIG.  4-20  BUTTERWORTH  SHAPED  DATA,  GAUSSIAN  PROBABILITY  DISTRIBUTION, 
1007.  DUTY  CYCLE 
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One  Slnewave  plus  Butterworth  Gaussian  Data  Test  Results 
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Figures  4-21  through  4-24  show  the  test  results  for  a  data  signal 
comprising  a  sinewave  plus  Butterworth  shaped  Gaussian  noise. 

The  tests  were  conducted  for: 

1)  various  values  of  the  ratio  of  the  data  cutoff  frequency 
to  the  sinewave  frequency,  fj/fj^ 

2)  two  values  in  db  of  the  ratio  of  the  sinewave  voltage  to 
the  noise  voltage,  As/Aj^  =  0  and  12  db 

3)  “  0.02  pulsewidth  and  100%  duty  cycle  pulsewidth. 

The  BWT  interpolation  filter'  was  used  throughout  these  tests. 

It  will  be  noticed  that  for  ratios  of  sinewave  to  noise  voltage 
of  12  db  the  test  results  are  similar  to  the  single  sinewave  test  results. 
At  ratios  of  zero  db  test  results  reflect  the  influence  of  the  Butterworth 
noise  in  that  the  minimum  rms  error  is  higher. 

The  100%  duty  cycle  interpolation  competes  favorably  with  the 
tiff  X  0.02  pulsewidth  case.  In  some  cases,  fj/fi  ®  2,  the  100%  duty 
cycle  pulsewidth  gives  the  best  results.  For  fj/fj^*!  the  =  0.02 

pulsewidth  sampling  is  better.  However,  the  differences  are  never  large. 

In  general,  as  the  sinewave  is  reduced  in  frequency  lower  rms 
error  is  achieved  for  a  fixed  sample  per  cycle  of  the  Gaussian  noise  cutoff 
frequency.  This  is  due  to  the  fact  that  the  sinewave  is  the  predominant 
data  and  as  its  frequency  is  reduced  the  samples  per  cycle  of  the  sinewave 
increases . 

Three  Sinewaves  plus  Butterworth  Gaussian  Data  Test  Results 

Figures  4-25  through  4-28  show  results  for  three  sinewaves  plus 
Butterworth  shaped  Gaussian  noise.  The  sinewaves  are  located  at  various 
points  in  the  data  spectrum:  fj/^l  “  ^d^^2  “  ^d/%  ~  2.  The  db 

ratio  of  the  sinewave  to  noise  voltage  is  variable  from  figure  to  figure. 
Thu  two  curves  in  each  figure  are  for  0.02  pulsewidth  and  100%  duty 

cycle  pulsewidth. 
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RMS  ERROR,  % 
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SAMPLES  PER  CYCLE,  fg/fj 


FIG.  4-23  BUTTERWORTH  PLUS  SINEWAVE  DATA,  =  0.02,  Ag/A  =  12  db 
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RMS  ERROR 


0  4  8  12  16  20 

SAMPLES  PER  CYCLE,  fg/fj 


FIG.  4-24  BUTTERWORTH  PLUS  SINEWAVE  DATA,  =  0.02,  Ag/Aj,j  =  0  db 
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RMS  ERROR 


FIG.  4-25  BUTTERWORTH  PLUS  THREE  SINEWAVES  DATA,  Ai /Aj^-A2/Aj^A3/Aj^0  db 
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RMS  ERROR 


FIG.  4-26  BUTTERWORTH  PLUS  THREE  SINEWAVES  DATA,  A,/a«=6  db  A  /Am=  0  db 
A3/AN  «=  12  db  N  >  2  N  . 
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RMS  ERROR 


FIG.  4-27 


BUTTERWORTH  PLUS  THREE  SINEWAVES  DATA,  Ai/Af^A2/Afj=A3/A  =12  db 
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RMS  ERROR 


FIG.  4-28  BUTTERWORTH  PLUS  THREE  SINEWAVES  DATA, 

Ai/Aj^A2/Afj»A3/Ajj=12  db;  i'i/fjj  =  2 
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Figures  4-25  and  4-27  are  nearly  the  same.  It  follows  that  even 
for  the  zero  db  case  the  noise  Is  relatively  unimportant  In  determining  the 
error.  Figure  4-28  displays  lower  rms  errors  for  a  given  sample  per  cycle 
of  the  data  cutoff  frequency  since  the  slnewave  frequencies  have  been 
lowered. 

Step  Type  Butterworth  Gaussian  Data 

The  test  results  and  the  spectrum  are  presented  in  Figure  4-29. 

The  results  show  that  the  low  frequency  data  is  predominant  in  determining 
the  aliasing  error.  As  the  ratio  of  A]^  to  A2  is  decreased  the  results  will 
approach  that  for  Butterworth  Gaussian  data. 

Bandpass  Data 

Figure  4-30  shows  the  results  for  bandpass  data  with  two  pairs 
of  upper  and  lower  cutoff  frequencies.  Values  of  tj^f^  "  0.02  and  100%  pulse- 
width  were  used.  The  BWT  Interpolation  filter  was  used.  Best  results  were 
obtained  when  the  low  frequency  data  signal  was  increased. 

Some  bandpass  Interpolation  tests  not  reported  were  tried  briefly, 
but  performance  was  poorer  than  for  lowpass  interpolation. 

4.3.3  Time  Function  Photographs  and  Spectral  Measurements 

The  following  data  is  included  to  illustrate  the  effects  observed 
in  performance  of  the  experimental  tests  and  typical  power  spectral  density 
measurements  taken  during  the  course  of  the  program. 

Time  Function  Photographs 

Figure  4-31  provides  a  comparison  of  the  data,  100%  duty  cycle 
pulsetrain,  and  the  BWT  Interpolation  filter  output  waveforms  for  equi- 
probable  amplitude  data.  The  delay  between  the  data  and  interpolation  fil¬ 
ter  output  is  accounted  for  by  the  analog  error  comparator. 

Figure  4-32  is  a  comparison  of  the  data  and  BWT  interpolation 
filter  output  waveforms  for  Gaussian  data  and  fractional  duty  cycle  pulse- 
train.  In  Figure  4-32  the  aliasing  spectrum  feedthrough  is  quite  evident 
in  the  second  photograph  from  the  top,  and  contributes  significantly  to  the 
error. 
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RMS  ERROR, 


04  8  12  16  20 

SAMPLES  PER  CYCLE,  fg/fj,  (fd=f2) 


FIG.  4-29  CONTINUOUS  WITH  STEP  BUTTERWORTH  SHAPED 
GAUSSIAN  DATA,  1007.  DUTY  CYCLE 
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RMS  ERROR 


FIG.  4-30  BANDPASS  GAUSSIAN  DATA  WITH  24  db/oct.  SLOPES 
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INTERPOLATION  FILTER  OUTPUT 


DATA 


DATA 

SAMPLED  DATA  OUTPUT 

J-7,  RMS  ERROR  =2.87. 


FIG.  4-31  PHOTOGRAPHS  OF  DATA,  BVfT  INTERPOLATION  FILTER  OUTPUT, 
AND  SAMPLER  OUTPUT.  EQUI -PROBABLE  AMPLITUDE 
DISTRIBUTION.  lOOX  DUTY  CYCLE 
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J.8,  fi/fd-2 


INTERPOLATION  FILTER  OUTPUT 
DATA 

fg/fj-  14,  RMS  ERROR  -  0.2% 


INTERPOLATION  FILTER  OUTPUT 
DATA 

fg/fd"  ^  error  -  9.4% 


INTERPOLATION  FILTER  OUTPUT 
DATA 

fs/fd"  3,  RMS  ERROR  -  2.6% 


FIG.  4-32  PHOTOGRAPHS  OF  DATA  AND  BWT  INTERPOLATION  FILTER  OUTPUT. 
GAUSSIAN  PROBABILITY  DISTRIBUTION  tifi-0.02. 
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J-8,  fi/fd-2 


J-7, 


J-7,  fi/fd-1 


INTERPOLATION  FILTER  OUTPUT 
DATA 

fg/fj-  14,  RMS  ERROR  -  0,847. 


INTERPOLATION  FILTER  OUTPUT 
DATA 

fg/fd-  2,  RMS  ERROR  -  4,257. 


INTERPOLATION  FILTER  OUTPUT 
DATA 

fg/fd-  8,  RMS  ERROR  -  3,757. 


FIG.  4-32  PHOTOGRAPHS  OF  DATA  AND  BWT  INTERPOLATION  FILTER  OUTPUT, 
GAUSSIAN  PROBABILITY  DISTRIBUTION  tifi-0.02. 
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Spectral  Measurements 


The  accompanying  spectral  measurements  have  been  recorded  using 
a  B&K  Audio  Frequency  Spectrum  Recorder  Type  3311  which  consists  of  AF 
Spectrometer  Type  2111  and  Level  Recorder  Type  2304.  It  Is  designed  for 
the  automatic  recording  of  spectrograms,  and  has  a  frequency  range  of 
50  cps  -  35  kcps.  The  bandpass  filters  are  either  full  octave  or  1/3 
octave  wide,  a  constant  percentage  of  the  center  frequency,  which 
produces  a  positive  6  db/oct.  change  In  slope  In  all  spectra  analyzed. 

This  slope  Is  readily  seen  In  Figure  4-33.  The  fact  that  the 
bandwidth  Is  broad  means  that  the  resulting  spectral  measurements  are 
lacking  In  detail.  This  Is  evident  in  Figure  4-33  for  the  equl-probable 
amplitude  noise  source  which  has  a  sin  (x)/x  spectrum  with  the  first  null 
at  4000  cps.  This  spectrum  is  essentially  flat  up  to  200  cps. 

The  spectra  of  Interest  In  the  single  channel  system  are  shown 
In  Figure  4-34.  That  portion  of  the  spectrum  which  contributes  to  aliasing 
Is  readily  seen  as  that  portion  above  500  cps  In  the  Sampled  Data  Spectrum. 

Figure  4-35  shows  comparable  spectra  In  the  single  channel  system 
for  the  equl-probable  amplitude  data  source. 

The  Slgmatron  noise  generator  provides  three  bandwidth  selec¬ 
tions  -  1000,  100,  and  10  cps.  When  the  data  bandwidth  was  200  cps  the 
1000  cps  bandwidth  was  used  and  the  resulting  cutoff  rate  was  that  due 
to  the  filter,  42  db/oct.  For  a  data  bandwidth  of  100  cps  the  Slgmatron 
noise  generator  bandwidth  was  set  at  100  cps  In  order  to  provide  sufficient 
signal  at  the  data  filter  output  and  the  cutoff  rate  of  the  noise  generator 
at  100  cps  is  approximately  6  db/oct.  Combining  this  cutoff  with  that  of 
the  filter  results  in  an  overall  rolloff  of  48  db/oct.  These  cutoff  rates 
were  measured  with  a  General  Radio  Company  Wave  Analyzer,  Type  736-A  and 
are  shown  in  Figure  4-36.  The  Slgmatron  noise  generator  output  for  100  cps 
bandwidth  setting  was  measured  as  well  as  the  filtered  output  using  the 
SKL  Model  308  filter.  When  a  data  bandwidth  of  50  cps  was  desired  the 
100  cps  Slgmatron  bandwidth  was  used  and  the  resulting  data  cutoff  rate 
was  again  42  db/oct.  This  problem  was  not  encountered  for  equl-probable 
amplitude  distribution  noise  since  the  3  db  noise  bandwidth  was  approximately 
ten  times  the  largest  data  bandwidth. 
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FIG.  4-34  SINGLE  CHANNEL  SYSTEM  SPECTRA, 
DATA:  GAUSSIAN  AMPLITUDE 
PROBABILITY  DISTRIBUTION 
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FIG.  4-36  SPECTRA  OF  BUTTERWORTH  SHAPED  DATA, 

GAUSSIAN  AMPLITUDE  PROBABILITY  DISTRIBUTION 
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FREQUENCY,  CPS 


4.4  MULTI-CHANNEL  SYSTEM  TESTS 


4.4.1  System  Description 

The  general  system  block  diagram  and  description  is  given  in 
Figure  5-1  Section  5.  A  brief  description  of  the  system  operation  is 
given  here. 

The  test  setup  is  similar  to  that  used  for  the  PAM-FM  tests 
described  in  Reference  1,  Volume  III.  The  same  subcarrier  equipment  was 
used  as  an  RF  carrier.  All  tests  were  performed  with  150  second  pulse- 
width  sampling  constituting  60Z  of  the  pulse  period  for  driving  the  inter¬ 
polation  filter.  The  analog  error  comparator  was  activated  only  by  the 
data  channel  being  measured  with  modulation  on  all  other  channels  removed 
to  permit  ready  determination  of  the  crosstalk  error.  No  pulse  detector 
was  used.  The  equivalent  function  was  performed  by  the  interpolation 
filter  since  the  wide  sample  pulsewidths  including  noise  were  fed  to  the 
lowpass  Interpolation  filter  which  discriminated  against  the  noise  by  per¬ 
forming  an  integrating  action  similar  to  that  performed  by  an  integrating 
pulse  detector.  The  position  of  the  decoimnutator  gate  within  the  pulse 
period  was  also  adjustable.  The  position  was  chosen  to  give  the  minimum 
system  residual  error  reading.  This  position  was  approximately  at  the 
center  of  the  pulse. 

4.4.2  Test  Results 


The  multi-channel  tests  served  to  establish  that  overall  system 
error  could  be  obtained  by  combining  Interpolation  error  determined  from 
single  channel  tests  with  the  other  error  sources  consisting  primarily  of 
interference  and  crosstalk.  Error  due  to  interference  and  crosstalk  were 
obtained  during  the  telemetry  system  study  reported  in  Reference  1.  This 
was  verified  for  the  P.'.M-FM  system. 

The  carrier  system  was  set  up  as  follows; 

1)  Receiver  bandwidth,  B  •  12  kcps 

2)  System  pulse  rate,  Fg  -  4  kcpa 

3)  M  -  B/Fg  -  3 

4)  Channel  sampling  frequency,  fg  =  500  cps 


87 


5)  Peak  deviation  -  +  4.5  kcps 

6)  No  premodulation  filtering 

7)  Lovpass  filter  following  discriminator  -  6  kcps 

8)  Data  -  Butterwrth,  Gaussian  noise 

9)  Interpolation  Filter  -  Butterworth 

Four  tests  were  conducted  with  this  setup.  Slope  of  the  data 
and  the  interpolation  filter  were  selected  and  tests  conducted  for  variable 
samples  per  cycle.  The  rms  error  was  measured  as  a  function  of  S/N,  the 
signal  to  noise  ratio,  measured  preceding  the  discriminator.  The  results 
of  these  tests  are  shown  in  Figures  4-37  through  4-40.  Calculated  results 
are  also  shown. 

The  calculated  results  were  obtained  by  adding  the  mean-square 
interpolation  error  for  the  single  channel  tests  to  the  mean-square  pulse 
error  with  the  residual  equipment  error  removed  from  the  values  of  S/kj^Fg^ 
obtained  from  Figure  11-6-17,  Reference  1,  Volume  II  reproduced  in  Figure 
4-41.  The  rms  error  was  obtained  by  taking  the  square  root  of  this  sum. 

The  agreement  obtained  is  very  good  in  most  cases.  This  method  allows  the 
overall  system  performance  to  be  calculated  for  any  of  the  time-division 
multiplexed  systems  under  a  wide  range  of  operating  conditions  and  signal- 
to-noise  ratios  by  using  the  single  channel  interpolation  error  in  con¬ 
junction  with  the  system  performance  determined  from  Reference  1  or  from 
other  data. 

4.4.3  Time  Function  Photographs  and  Spectral  Measurements 

Time  Function  Photographs 

Figure  4-42  shows  multiple  sweep  oscilloscope  waveforms  at  the 
8  channel  multiplexer  output.  The  top  photograph  shows  several  frames 
with  all  channels  modulated;  the  next  photograph  has  an  expanded  sweep 
rate  to  show  detail  as  does  the  bottom  photograph  where  only  one  channel 
is  modulated. 

Figure  4-43  shows  the  decommutator  input  and  output  waveforms 
with  and  without  adjacent  channel  modulation.  It  can  be  seen  that  the  out¬ 
put  waveform  is  a  gated  portion  of  the  input  waveform. 

Spectral  Measurements 

Spectra  were  observed  at  various  points  in  the  system  and  were 
found  to  appear  similar  to  those  shown  for  the  single  channel  tests  at  the 
corresponding  points.  Carrier  and  multiplex  output  spectra  were  checked  and 
are  as  reported  in  Reference  1. 
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FIG.  4-37  RMS  ERROR  VS.  S/kjFg^.  INTERPOLATION  FILTER:  BUTTERWORTH, 
24  db/oct.  DATA:  BUTTERWORTH,  J  -  4  AND  5 
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RMS  ERROR 


S/kjFg^,  db 


FIG.  4-38  RMS  ERROR  VS.  S/kjF,^.  INTERPOLATION  FILTER:  BUTTERWORTH 
24  db/oct.  DATA:  BUTTERWORTH,  J  -  7  AND  8 
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S/k^F,^ 


FIG.  4-39  RMS  ERROR  VS.  S/kiFg^.  INTERPOLATION  FILTER:  BUTTERVIORTH, 
36  db/oct.  DATA:  BUTTERVWRTH,  J  -  4  AND  5 
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RMS  ERROR 


FIG.  4-AO  RMS  ERROR  VS.  S/k,F  INTERPOLATION  FILTER:  BUTTERWORIH, 
36  db/oct.  DATA;  iuTTERMORTH,  J  «  7  AND  8 
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RMS  ERROR 


FIG.  4-41  PAM-FM,  RMS  ERROR  VS.  S/k,F  RMS  DEVIATION, 

fD/h  -  0-A. 
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%  FRAME 

1  CHANNEL  MODULATED 


FIG.  4-42  PHOTOGRAPHS  OF  MULTIPLEXER  OUTPUT  WAVEFORM 
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SECTION  5 


TEST  EQUIPMENT 


5 . 1  INTRODUCTION 

The  test  equipment  designed  to  measure  experimentally  the 
aliasing  error  In  telemetry  systems  Is  described  In  this  section.  The 
general  test  setup  will  be  presented  In  Section  5.2.  This  will  Include 
the  single  channel  and  multi-channel  descriptions.  Section  5.3  will 
discuss  the  equipment  making  up  the  telemetry  system.  The  analog  error 
comparator  Is  described  In  Section  5.4,  and  the  accuracy  with  \dilch 
measurements  were  obtained  Is  evaluated  In  Section  5.5. 

Many  of  the  basic  elements  of  the  system,  the  analog  error 
comparator,  and  the  evaluation  of  the  accuracy  of  measurements  have 
been  previously  discussed  In  reference  1,  Volume  III.  Reference  will 
be  made  to  this  volume  frequently.  Therefore,  It  will  be  referred  to 
simply  as  "reference  1". 

The  majority  of  the  equipment  used  for  this  contract  was  trans¬ 
ferred  from  the  Telemetry  System  Study,  Contract  No.  DA-36-039-SC-73182 , 
DA  Project  No.  3-16-00-300.  Two  additional  pieces  of  equipment  were  de¬ 
signed  for  this  project;  1)  an  8-channel  electronic  multiplexer,  and  2) 
an  electronic  decommutator.  These  are  described  In  Section  5.3. 


5,2  GENERAL  TEST  SETUP 

The  setup  for  measuring  errors  due  to  aliasing  and  Interpolation 
for  a  single -channel  or  time-division  multiplex  sampled  data  system  is 
shown  in  Figure  5-1.  Basically  the  setup  provides  an  assembly  of  a  typical 
telemetry  system  Including  a  signal  source  which  simulates  the  telemetry 
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FIG.  5-1  BLOCK  DIAGRAM  -  GENER^  TEST  SETUP 


data,  sources  of  Interference,  monitoring  equipment,  and  an  analog  error 
comparator . 

The  test  setup  provided  for  adjustment  of  significant  parameters 
in  order  that  a  thorough  evaluation  could  be  accomplished.  Signal  simula¬ 
tors  are  multiplexed  to  provide  the  modulating  wave  to  the  telemetry  sys¬ 
tem.  The  duty  cycle  of  the  multiplexer  and  decommutator  output  can  be 
varied.  The  error  measurement  made  by  the  amilog  error  comparator  Is 
accomplished  by  comparing  the  signal  simulator  output  with  the  interpola¬ 
tion  filter  output.  The  data  signal  is  taken  to  be  the  data  filter  output 
since  the  filter  was  utilized  to  simulate  typical  data  spectra  by  shaping 
broadband  noise  sources. 

In  order  to  evaluate  the  effects  of  crosstalk  three  signals  can 
be  multiplexed  at  the  input  to  the  communication  link.  Separate  data  sig¬ 
nals  modulate  the  measured  channel  and  the  multiplexed  channel  preceding 
the  measured  channel.  The  third  signal  is  a  quasi-random  PAM  pulsetrain 
filling  up  the  6  remaining  channel  positions.  Separate  control  of  the 
channel  preceding  Che  measured  channel  allows  determination  of  crosstalk 
effects . 


The  single  channel  sampled  data  system  is  obtained  by  omitting 
the  multiplexer  and  the  carrier  equipment.  The  experimental  setup  con¬ 
sisted  of  the  following  elements  in  tandem:  data  source,  data  filter, 
decommutaCor ,  interpolation  filter,  and  amplifier.  The  error  measure¬ 
ments  were  made  as  before  by  comparing  the  data  filter  output  with  the 
amplified  interpolation  filter  output. 


5.3  SYSTEM  EQUIPMENT 

The  majority  of  the  equipment  used  for  the  aliasing  project  has 
been  discussed  in  reference  1.  Rather  than  repeating  descriptions  of 
b-sic  elements  they  will  be  itemized  and  the  reader  may  obtain  a  detailed 
description  there. 

5.3.1  Signal  Simulators 

1)  Random  Noise  Generator,  Model  Ot  1390-A,  General  Radio  Co. 
(Gaussian  Amplitude  Probability  Distribution). 

2)  Random  Signal  Generator,  Sigmatron,  Model  GA-1000,  Inter¬ 
continental  Dynamics  Corporation,  (Gaussian  Amplitude 
Probability  Distribution). 
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3)  Slnewave  Oscillator,  Model  HP-200CD,  Hewlett-Packard  Co. 

4)  Random  Noise  Generator,  Modified  Navcor  Equipment,  (Equal 
Amplitude  Probability  Distribution). 


5.3.2  Data  Filter 

1)  Variable  Filter,  Model  308,  Spencer  Kennedy  Laboratories, 

Inc . 

The  frequency  response  for  this  filter  is  shown  in  Figures 
5-2  and  5-3  for  one  and  two  sections  in  cascade  respectively. 

A  frequency  response  curve  is  shown  for  each  of  the  data 
bandwldths  considered  in  the  experimental  test  program.  The 
cutoff  rates  were  determined  at  the  30  db  attenuation  point. 
One  section  of  the  filter  gives  a  cutoff  rate  of  24  db/oct., 
while  two  sections  in  cascade  provide  a  42  db/oct.  cutoff 
rate  (the  specifications  quote  48  db/oct.). 


5.3  3  Multiplexer 


General  Description 

The  multiplexer  provides  a  PAM  wave  of  either  507.  or  1007.  duty 
cycle.  Between  pulses  the  507.  duty  cycle  wave  returns  to  zero  or  to  any 
desired  pedestal  level. 

Three  Independent  signal  sources  may  be  multiplexed  to  provide 
eight  PAM  pulses  per  frame.  Channels  1  and  2  are  independent,  while  chan¬ 
nels  3  through  8  are  strapped  at  the  input.  Channel  2  is  used  for  test. 

An  external  sine  wave  clock  controls  the  PAM  pulse  rate.  The 
multiplexer  will  operate  at  rates  from  500  to  5000  pulses  per  second.  Full 
scale  PAM  output  amplitude  is  ±  2.5  volts.  A  frame  sync  pulse  is  provided 
for  the  decommutator. 

The  multiplexer  is  shown  in  block  form  on  Fig.  5-4. 
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FIG.  5-2  FILTER  RESPONSE  -  SKL  MODEL  308,  24  DB/OCT. 


100 


500 


FIG.  5-3  FILTER  RESPONSE  -  SKL  MODEL  308,  42  DB/OCT. 
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FIG,  5-4  MULTIPLEXER  BLOCK  DIAGRAM,  AL-1 


Circuit  Description 

Figures  5-5  through  5-11  show  circuit  details  schematically. 

The  Input  amplifier.  Fig.  5-5,  is  an  operational  type.  Two  in¬ 
puts  are  provided  so  that  It  may  be  connected  to  Invert  or  not  Invert 
the  signal.  In  the  multiplexer  no  signal  Inversion  takes  place  between 
Input  and  output. 

The  coninutator  switch.  Fig.  5-6,  Is  of  the  series  type.  The 
gating  transistors  are  driven  to  saturation  so  the  effective  series  Im¬ 
pedance  Is  low.  The  ratio  of  the  off  to  the  on  Impedance  Is  greater 
than  105.  This  type  of  switch  causes  a  bias  current  in  the  input  cir¬ 
cuit  so  it  must  be  driven  by  a  low  Impedance  source.  The  operational 
amplifier  Is  satisfactory  for  this. 

The  outputs  of  the  three  commutator  switches  are  combined  at 
the  input  of  another  operational  amplifier.  The  multiplexed  signals  are 
then  sent  to  the  gate  and  holding  circuit.  Fig.  5-7.  If  a  PAM  wave  of 
1007.  duty  cycle  is  required,  the  signal  is  sampled  for  approximately 
20  >A  sec  during  each  channel  period  and  then  held  until  the  next  sample. 

A  non-return  to  zero  PAM  wave  is  generated.  If  a  507.  duty  cycle  wave  is 
required,  the  sample  time  is  extended  to  one-half  of  the  channel  period 
but  no  holding  action  is  performed.  During  the  time  between  samples  the 
PAM  wave  returns  to  an  adjustable  pedestal  value.  The  gate  and  holding 
circuit  contains  a  series  gate  and  an  output  amplifier. 

The  remainder  of  the  multiplexer  generates  the  timing  waves 
that  control  the  commutator  gates  and  the  ouput  gate. 

A  trigger  generator.  Fig.  5-8,  provides  sharp  pulses  to  drive 
the  counter.  An  Inductive  kickback  pulse  generator  is  used.  One  out¬ 
put  pulse  is  generated  per  cycle  of  input  from  the  sine  wave  clock. 

A  binary  counter,  consisting  of  three  flip  flops  (Fig.  5-9), 
provides  three  waves  which  are  combined  in  the  multiplexer  gate  genera¬ 
tor,  Fig,  5-10,  to  generate  the  three  sequential  timing  waves  to  control 
the  commutator  switches.  The  control  wave  for  Channel  1  is  also  ampli¬ 
fied  and  used  as  a  frame  sync  pulse  for  deconmutation.  The  amplifier 
for  this  is  shown  on  Fig.  5-8. 

An  additional  timing  wave  is  generated  by  the  sample  gate 
generator,  Fig.  5-11.  The  sampling  gate  is  delayed  approximately  50/1  sec 
from  the  beginning  of  each  channel  period  to  ensure  that  the  commutator 
switching  transients  have  died  away  before  sampling  of  the  combined  wave 
takes  place. 
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FIG.  5-5  AMPLIFIER 
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FIG.  5-6  COMMUTATOR  SWITCH,  AL-4 
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FIG.  5-8  TRIGGER  GENERATC® 
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FIG.  5-9  FLIP-FLOP,  AL- 


FIG.  5-11  SAMPLE  GATE  GENERATOR 


Three  control  functions  are  provided  on  the  multiplexer:  duty 
cycle,  pedestal  level,  and  sample  pulse  duration. 


Operating  Adjustments 

1.  Set  the  clock  generator  to  a  frequency  corresponding  to  the 
required  channel  rate.  The  clock  wave  amplitude  must  be  at  least  10 
volts  rms. 

2.  Adjust  the  input  signal  levels  so  that  the  amplitude  of  the 
output  PAM  wave  is  +  2.5  volts  full  scale  on  each  channel. 

3.  Set  the  duty  cycle  selector  switch  to  the  50%  position.  Ad¬ 
just  the  coarse  and  vernier  gate  width  controls  for  the  required  gate 
width.  If  100%  duty  cycle  PAM  is  required  adjust  the  width  controls  for 
a  20  ya  sec  gate  and  then  switch  to  100%  duty  cycle. 

4.  If  50%  duty  cycle  PAM  is  required,  adjust  the  pedestal  level 
control  for  the  desired  pedestal  level. 
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5.3.4  Conanunlcatlon  Link 


The  coninunlcatlon  link  shown  in  Figure  5-1  consists  of  the 

following: 

1)  Amplifier 

2)  Subcarrier  Oscillator,  Bendix  TGO-601 

3)  Bandpass  Filter  and  Mixing  Amplifier 

4)  Subcarrier  Discriminator,  EMR  97D 

5)  Adjustable  Lowpass  Filter,  EMR  95D 
See  Reference  1  for  equipment  descriptions. 

5.3.5  Decommutator 


General  Description 

The  decommutator  provides  a  means  of  separating  the  desired 
channel  from  the  reamlnder  of  the  multiplexed  PAM  wave.  Either  of  two 
output  gates  may  be  used.  One  provides  a  variable  aperture  sample  with 
return  to  zero  and  the  other  a  short  sample  which  is  held. 

Two  delay  circuits  are  provided.  One  delays  the  input  PAM 
wave  so  that  the  analog  error  comparator  may  be  used  to  measure  the  error 
between  the  input  and  output  of  the  deconmutator .  The  other  delays  the 
frame  sync  pulse  from  the  multiplexer  and  can  be  adjusted  so  that  the  de- 
commutator  gate  operates  at  the  proper  time. 

A  frequency  divider  is  provided  to  select  only  every  second, 
third,  or  fourth  frame  sync  pulse  in  case  the  analog  error  comparator  is 
to  be  operated  synchronously  with  the  multiplexer  at  high  pulse  rates. 

The  divider  is  not  used  when  the  analog  error  comparator  samples  asynch¬ 
ronously. 


Controls  are  provided  to  adjust  the  input  bias,  the  output  bias, 
the  sampling  pulse  delay,  and  the  sampling  pulsewldth. 

Tlie  decommutator  block  diagram  is  shown  in  Figure  5-12. 
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FIG.  5-12  DECOMMUTATOR  BLOCK  DIAGRAM,  AL-2 


Circuit  Description 


Figures  5-13  througn  5-18  show  circuit  details  schematically. 

The  delay  network.  Fig.  5-13,  is  an  11-pole  maximally  flat  time 
delay  type.  It  has  a  1  db  bandwidth  of  500  cps, which  is  several  times  the 
maximum  data  bandwidth,  and  a  time  delay  of  approximately  600 /i.  sec  re¬ 
quired  for  proper  operation  of  the  error  comparator.  In  the  multichannel 
tests  the  delay  network  is  removed  because  the  carrier  system  supplies  this 
delay. 


An  amplifier.  Fig.  5-5,  drives  the  decomnutator  gate.  A  switch 
permits  polarity  inversion  of  the  PAM  wave  if  it  is  required.  This  ampli¬ 
fier  is  the  same  as  those  used  in  the  multiplexer. 

The  first  type  of  decommutator  gate,  Fig.  5-14,  provides  variable 
aperture  gating  of  the  PAM  wave. 

It  contains  an  amplifier  capable  of  driving  a  2k  ohm  LC  inter¬ 
polation  filter.  The  gating  circuit  is  similar  to  that  in  the  commutator 
switch  in  the  multiplexer.  The  gate  output  returns  to  zero  or  to  a  se¬ 
lected  bias  value  after  sampling.  The  second  type  of  gate,  Fig.  5-15,  is 
similar  except  that  it  provides  a  holding  circuit  which  maintains  the  out¬ 
put  at  the  previous  sample  level  until  the  next  sample  begins. 

The  remainder  of  the  decommutator  circuitry  provides  a  sample 
gate  pulse  to  the  decommutator  gate  at  the  proper  time. 

The  trigger  generator.  Fig.  5-16,  generates  sharp  pulses  from 
either  a  sine  wave  clock  or  from  the  frame  sync  pulse  from  the  multiplexer. 

The  clock  pulses  are  delayed  by  two  adjustable  one  shot  multi¬ 
vibrators,  Fig.  5-17.  A  third  one  shot  generates  the  sample  gate  pulse. 
Coarse  and  vernier  controls  are  provided  for  each  function. 

9 

The  frequency  divider.  Fig.  5-18,  is  a  synchronized  multi¬ 
vibrator.  Pulses  from  the  trigger  generator  are  injected  to  lock  the 
period  of  the  multivibrator  wave  at  an  integral  multiple  of  the  period  of 
the  clock  pulse.  An  injection  sensitivity  control  and  coarse  and  vernier 
frequency  controls  are  provided. 


FIG.  5-15  DECOMMUTATOR  GATE  AND  HOLDING  CIRCUIT 
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FIG.  5-16  TRIGGER  GENERATOR,  AL-13 
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FIG.  5-18  FREQUENCY  DIVIDER,  AL-15 


Operating  Adlugtments 


1.  Adjust  the  sine  wave  clock  (if  frame  sync  from  the  multi¬ 
plexer  is  not  used)  to  the  required  sample  frequency  between  100  and  2000 
cps.  The  amplitude  must  be  at  least  10  volts  rms. 

2.  Observe  the  decommutator  output  with  an  oscilloscope.  Ad¬ 
just  the  amplifier  bias  and  the  PAM  Input  level  so  that  the  decommutator 
output  is  +  2.5  volts  full  scale. 

3.  Adjust  the  output  bias  control  so  that  the  output  level 
between  samples  is  zero. 

4.  Insert  AL-11  In  the  decommutator  gate  socket.  Adjust  the 
sample  gate  width  controls  for  the  required  aperture. 

5.  Adjust  the  sample  delay  controls  so  the  gate  Is  positioned 
properly  In  the  frame  for  correct  decommutatlon.  If  the  holding  circuit 
Is  desired,  replace  AL-11  with  AL-12. 

6.  If  synchronous  operation  of  the  analog  error  comparator  Is 
required,  adjust  the  divider  sensitivity  and  frequency  controls  until  a 
stable  wave  is  obtained  at  the  desired  frequency. 
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5.3.6  Interpolation  Filters 


Four  Interpolation  filters  were  Investigated.  These  are: 

1)  A  200  cps  cutoff  filter  with  a  trap  at  400  cps  with  the 
frequency  response  whown  In  Figure  5->19.  This  Is  a  6  element  Butter- 
worth  filter  with  an  added  capacitor  to  form  a  trap  with  the  center  In 
doctor.  The  resulting  cutoff  rate  Is  approximately  60  db/oct.  at  the 
30  db  attenuation  point. 

2)  A  Gaussian  interpolation  filter  with  a  final  cutoff  rate 
of  approximately  24  db/oct.  The  frequency  response  Is  shown  in  Figure 
5-20. 


3)  A  4-pole  Butterworth  filter  with  the  frequency  response 
shown  In  Figure  5-21. 

4)  A  6-pole  Butterworth  filter  with  the  frequency  response 
shown  In  Figure  5-22. 
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FIG.  5-19  FILTER  RESPONSE  -  BWT,  60  DB/OCT. 
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RESISTANCE  IN  OHMS 


FIG. 
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5-20  FILTER  RESPONSE  -  EMR  95D,  GAUSSIAN  CHARACTERISTIC,  24  DB/OCT. 
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FIG.  5-21  FILTER  RESPONSE  -  4-POLE  BUTTERWORTH  FILTER,  24  DB/OCT. 
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FREQUENCY,  CPS 


FIG.  5-22  FILTF.R  RESPONSE  -  6-POLE  BUTTERWORTH  FILTER,  36  DB/OCT, 
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FREQUENCY,  CPS 


5.4  ANALOG  ERROR  COMPARATOR 


Figure  5-23  shows  the  comparator  in  block  form.  A  detailed 
description  is  given  in  reference  1,  Vol.  III.  It  is  basically  a  device 
for  sampling  two  analog  signals  and  comparing  their  amplitudes.  The 
sampling  is  done  serially.  The  amplitude  differences  are  Indicated  In 
terms  of  both  average  and  rms  value. 

The  following  sequence  of  operations  takes  place  in  the  analog 
error  comparator: 

A.  The  sampling  switch  connects  the  follow  and  hold  to  the 
reference  analog  signal. 

B.  A  hold  pulse  causes  the  follow  and  hold  unit  to  hold  its 
output  signal  constant  for  200 yn.sec.  (sample  A). 

C.  A  convert  pulse  causes  the  Adage  A/D  converter  to  convert 
sample  A  into  digital  form  (parallel  binary). 

D.  A  store  pulse  causes  the  A/D  converter  to  put  this  binary 
number  into  a  storage  register. 

E.  The  sampling  switch  connects  the  follow  and  hold  unit  to 
the  output  of  the  telemetering  system. 

F.  A  second  hold  pulse  causes  the  follow  and  hold  unit  to 
hold  its  output  signal  constant  for  200  /X  sec.  (sample  B) . 

G.  A  second  convert  pulse  causes  the  A/D  converter  to  sub¬ 
tract  sample  B  from  the  analog  of  Sample  A,  convert  the 
difference  to  digital  form,  and  present  it  as  a  parallel 
binary  number  at  the  output  terminals. 

H.  A  dump  and  store  pulse  clears  the  storage  regi.'ster  in  the 
D/A  converter  and  enters  the  binary  number  from  the  A/D 
converter . 

J.  A  dump  pulse  clears  the  storage  register  in  the  A/D 
converter . 

This  completes  the  sequence  of  operation  of  the  analog  error 
comparator.  It  is  capable  of  making  up  to  2000  comparisons  per  second. 
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The  time  between  the  A  and  B  samples  is  determined  by  a  sine 
wave  clock.  Ihe  clock  frequency  Is  adjusted  so  that  the  period  of  one 
clock  cycle  equals  the  time  delay  in  the  telemetering  system. 


5.5  ACCURACY  OF  MEASUREMENTS 

Overall  accuracy  and  limitations  are  discussed  in  reference  1 
Volume  III. 

The  overall  measurement  accuracy  for  the  Aliasing  tests  are 
considered  to  be  within  the  limits  shown  below: 

Measurenent  Limiting  Accuracy 

Rms  error  57,  (0.05%  resolution) 

S/N  1  db 


Spectrum  amplitude 


3  db 


SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions  and  recommendations  derived  from  this  study 
are  summarized  below: 

6.1  CONCLUSIONS 

6.1.1  Determination  of  Typical  Data  Spectra 

1)  Satisfactory  spectral  models  have  been  derived  for  rep¬ 
resenting  typical  spectral  data  characteristics. 

2)  Practical  application  of  the  spectral  models  requires 
additional  information  on  the  part  of  the  system  designer.  This  informa¬ 
tion  consists  of  an  understanding  of  tie  basic  mechanisms  generating  the 
data  spectra  to  permit  specification  of  the  bandwidth  and  spectral  shapes 
of  the  data  spectra  by  the  designer. 

6.1.2  Analysis 

1)  The  analysis  and  numerical  computations  of  selected 
examples  agreed  closely  with  experimental  data. 

2)  Additional  numerical  computations  to  predict  aliasing  and 
interpolation  distortion  errors  should  be  undertaken.  These  computations 
should  be  based  on  the  spectral  models  derived  in  this  study  and  should 
take  into  account  a  broad  class  of  interpolation  filters  to  provide  a 
number  of  different  filter  types  to  instrumentation  engineers  for  re¬ 
constructing  data. 
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3)  An  analytical  Investigation  supplemented  by  Implementation 
of  a  practical  technique  for  obtaining  "allas>free"  sampling  should  be 
undertaken.  This  Investigation  should  be  supplemented  by  experimental 
test  and  evaluation. 

6.1.3  Laboratory  Measurements 

1)  The  experimental  test  program  conducted  In  this  study  has 
provided  proven  operating  parameters  for  a  time-division  telemetry  system 
which  guarantees  a  minimum  standard  of  performance. 

2)  Techniques  for  improving  this  performance  should  be  in¬ 
vestigated  experimentally,  especially  In  determining  the  error  mechanisms 
which  limit  the  ultimate  accuracy  of  sampled  data  systems.  The  range  of 
accuracies  which  should  be  investigated  lie  in  the  region  of  .05  to  .25% 
rms . 

6.1.4  Operating  Parameters 

The  most  salient  factor  In  controlling  Interpolation  error  Is  the 
spectral  characteristic  of  the  signal  to  be  sampled.  This,  of  course,  was 
well  known  and  is  borne  out  by  the  results  of  this  study.  This  means  that 
to  optimize  the  system  parameters  the  data  spectrum  must  be  known  and  de¬ 
sign  Information  such  as  given  in  Section  4  must  be  utilized  for  the  speci¬ 
fic  data  spectrum  being  sampled.  However,  when  the  data  spectrum  Is  not 
accurately  known  many  specific  recommendations  can  be  made  to  aid  the  sys¬ 
tem  design  based  upon  the  results  of  this  study. 

A  very  important  result  of  this  work  is  the  determination  of 
which  parameters  are  sensitive  In  controlling  aliasing  and  interpolation 
errors,  and  which  are  insensitive.  The  fact  that  performance  is  insensitive 
to  some  parameters  makes  the  design  information  in  Section  3  and  4  more 
universally  applicable. 

The  following  conclusions  relating  to  operating  parameters  can 
be  enumerated: 

1)  The  aliasing  error  is  primarily  dependent  upon  the  data 
spectrum  and  upon  fs/fdCthe  samples  per  cycle  of  data).  If  fs/fd  is  to 
minimized,  minimizing  the  required  power  and  bandwidth,  the  data  spectrum 
must  be  accurately  known , 

2)  The  amplitude  probability  distribution  of  the  data  does  not 
affect  the  aliasing  error  (or  distortion  error  wliich  is  due  to  amplitude 
and  time  delay  variation  with  frequency). 
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3)  It  is  the  high  frequency  region  of  the  data  spectrum  which 
is  of  most  importance  (unless  this  region  is  very  low  in  relative  energy) 
which  often  allows  simplifying  the  spectral  type  f'>r  ordinary  low-pass 
interpolation.  For  continuous  data  spectra  the  high  frequency  skirt  or 
rolloff  characteristic  is  very  important. 

4)  The  Interpolation  error  is  relatively  insensitive  to  the 

interpolation  filter  characteristic  for  the  low-pass  types  examined  with 
selectivity  from  24  to  60  db/oct.  The  most  selective  filter  performs 
best  at  low  samples  per  cycle,  2  to  3  or  4,  while  the  36  db/oct.  Butter- 
worth  type  is  about  optimum  for  higher  fg/fd  reduced  distortion 

error.  The  Gaussian  characteristic  is  inferior  to  the  Butterworth  for 
interpolation.  The  24  db/oct.  Butterworth  filter  performed  well  enough 
for  many  applications  and  would  be  useful  for  low  frequency  interpolation 
where  a  relatively  simple  R-C  mechanization  is  desired. 

5)  The  interpolation  error  is  insensitive,  to  the  duty  cycle 
of  the  interpolated  pulsetrain,  even  from  short  pulses  to  100%  duty  cycle 
for  the  same  interpolation  filters.  The  100%  duty  cycle  signal  supplies 
much  more  energy  to  the  interpolation  filter  minimizing  the  effects  of 
system  noise.  It  also  minimizes  the  modulation  effects  of  possible  pulse 
jitter.  This  waveform  will  therefore  improve  the  performance  of  most 
practical  systems  by  reducing  these  noise  sources. 

6)  In  most  applications,  especially  when  the  interpolation 
error  must  be  1%  or  better,  a  minimum  fg/fd  of  3  or  higher  will  be  re¬ 
quired.  For  fg/fd  ^  3  the  interpolation  filter  cutoff,  f^,  should  be  set 
for  fi/fd  of  about  2  which  minimizes  the  distortion  error  due  to  ampli¬ 
tude  and  phase  effects. 


6 . 2  RECOMMENDATIONS 

The  following  recommendations  are  based  on  the  conclusions 
cited  in  Section  6.1. 

1)  Investigation  of  the  nature  of  the  data  sources  should  be 
made  from  the  viewpoint  of  practical  instrximentation  applications.  Such 
an  investigation  would  consider  typical  data  sources  and  methods  of  es¬ 
timating  the  proper  bandwidth  and  spectral  shape  to  assign  to  the  given 
data  source.  Consideration  of  data  sources  derived  from  deterministic 
as  well  as  statistical  processes  and  the  common  limitations  imposed  on 
the  processes  by  physical  considerations  would  lead  to  more  realistic 
estimates  of  the  spectrum  of  the  data  sources.  A  simple  example  of  such 
estimation  techniques  would  be  to  take  velocity,  acceleration  and  posi¬ 
tion  limiting  into  account  in  the  derivation  of  the  spectrum  generated 
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by  a  physical  variable  such  as  temperature,  pressure,  shaft  positions, 
and  so  forth.  The  effects  of  data  spectral  filtering  by  pre-sampling 
filtering  could  then  be  more  realistically  evaluated  and  applied. 

2)  Computation  of  a  family  of  parametric  curves  presenting 
rms  error  versus  sampling  rate  for  different  spectral  data  models  and 
filter  types  should  be  undertaken.  This  data  should  be  based  on  selec¬ 
ted  sampling  and  reconstruction  techniques  representing  the  best  methods 
for  practical  Instrumentation  at  the  present  time. 

3)  An  experimental  investigation  of  techniques  to  obtain 
"alias-free"  sampling  should  be  undertaken. 

4)  Techniques  for  performing  more  accurate  data  sampling, 
processing  of  the  sampled  data  value  into  a  suitable  carrier  waveform, 
and  recovery  of  the  processed  data  sample  before  the  interpolation 
process  is  performed  should  be  undertaken  to  extend  the  accuracy  of  data 
recovery  into  the  accuracy  region  of  ,05  to  .25%  rms. 
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APPENDIX  1 


PAM  MULTIPLEX  SYSTEM  ANALYSIS 


Al-1  SUMMARY  OF  RESULTS 

This  Appendix  includes  the  following  results  applicable  to  the 
design  of  a  PAM  multiplex  system: 

a)  A  simplified  system  model  Including  the  effects  of  finite 
bandwidth,  finite  sampling  width,  aliasing,  crosstalk, 
and  noise  developed  for  an  N-channel  system. 

b)  A  general  expression  for  total  mean  square  error  due  to 
all  the  effects  In  a)  Is  obtained. 

c)  An  error  minimization  procedure  Is  carried  out  for  the 
special  case  In  which  all  channels  are  Identical  and  a 
signal  power  constraint  Is  imposed. 

d)  Expressions  for  optimum  transfer  functions,  both  realizable 
and  non-real Izable,  are  obtained  for  the  data  filter  and 
the  Interpolation  filter.  In  general,  these  transfer 
functions  are  bandllmlted  and  can  only  be  approximated  In 
the  practical  case. 

e)  Optimum  sampling  rate  is  specified  in  terms  of  the  filter 
transfer  functions. 

f)  Optimum  sampling  delay  Is  specified  In  terms  of  the  video 
filters  and  the  design  of  these  filters  to  reduce  crosstalk 
Is  mentioned. 
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Al-2  SYSTEM  MODEL  AND  NOTATION 


The  basic  system  model  v^lch  was  assumed  for  purposes  of  analy¬ 
sis  is  shown  in  Figure  Al-1.  This  model  Includes  N  input  channels  and 
one  output  channel. 


Xj.(t) 


FROM  TRANS-  VIDEO 

MISSION  CHANNEL  FILTER 


SAMPLING 

WINDOW 


SAMPLING  INTERPOLA-  1  OUT- 
TION  FILTER  PUT 


FIG.  Al-1  BLOCK  DIAGRAM  OF  N  CHANNEL  PAM  SYSTEM 


Tlie  various  transfer  functions  associated  with  one  channel  of 

this  model  are  discussed  in  detail  below.  Although  corresponding  transfer 
functions  for  the  other  channels  would  frequently  be  the  same,  this  con¬ 
dition  is  not  required  in  the  present  analysis.  It  may  also  be  noted  that 
modulation  to  an  RF  signal  and  back  again  is  not  considered. 
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H(](s)  >  The  data  filter  shapes  the  signal  spectrum  to  an 

optimum  form.  This  transfer  function  Is  assumed 
to  be  capable  of  adjustment. 

H^(s)  •  The  sampling  window  makes  It  possible  to  consider 

a  finite  sample  width.  Typically,  this  would  be 
the  transfer  function  for  a  finite-memory  inte¬ 
grator  In  which  the  Integration  time  Is  the  desired 
sample  width,  Tg.  Thus,  for  this  case 


Hw(8) 


1,  - 


fi(t) 


If  Instantaneous  samples  are  taken,  then 

Ms)  -  1 

The  actual  sampling  operation  is  represented  by  mul- 
tlpling  the  channel  signal  by  a  periodic  sequence 
of  Impulses  and  holding  the  product  for  a  specified 
time  interval.  For  the  nth  channel  the  sampling 
function  becomes 


fl(t)  -  X  ^  +  nt^-  kTj^)  ,  n.-  1,2, 


.,N 


k  »-*o 


%is) 


Where  Tj  is  the  sampling  interval,  per  channel,  and 
t]^  Is  the  Interval  between  sampling  Instants  In 
adjacent  channels. 

The  hold  network  would  usually  generate  rectangular 
pulses  from  the  Impulses.  If  the  pulse  width  Is 
tj^  (largest  possible  value)  then 

%(»)  “  1  - 


KyCs)  =  The  video  filter  smooths  the  boxcar  sequence  coming 
from  the  sunming  circuit  in  order  to  reduce  the 
channel  bandwidth.  The  minimum  bandwidth  is  deter¬ 
mined  by  crosstalk  considerations. 

Hy(s)  “  The  video  filter  at  the  receiver  input  is  similar  to 
the  one  above  and  is  used  to  reduce  the  noise  into 
the  sampling  circuits. 
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H^(s)  *  This  sampling  window  again  makes  It  possible  to  con¬ 
sider  a  finite  sample  width,  if  desired.  Its  form 
would  be  similar  to  H^(s). 

£2(t)  -  The  receiver  sampling  function  is  another  periodic 
sequence  of  impulses  given  by 

oO 

f2(c)  -  ^  (t  -  +  ntj  -  KTj)  ,  n  -  1,2 . N 

k  «  -*0 

Where  'C  is  the  sampling  delay  needed  to  compensate 
for  the  delay  in  the  transmission  channel. 

H^(s)  ■  The  interpolation  filter  converts  the  reconstructed 
samples  into  the  output  signal.  This  transfer 
function  is  assumed  to  be  capable  of  adjustment. 

In  addition  to  the  system  components  discussed  above,  noise  will 
be  Introduced  into  the  system  in  the  transmission  channel.  This  noise 
will  be  assumed  to  be  additive. 


Al-3  MODEL  SIMPLIFICATION 

For  purposes  of  notational  simplicity  it  is  desirable  to  represent 
any  one  channel  of  the  system  in  the  more  compact  form  shown  in  Figure 
Al-2.  The  transfer  functions  in  this  model  may  be  expressed  in  terms  of 
those 


FIG.  Al-2  SIMPLIFIED  BLOCK  DIAOtAM  FOR  A  SINGLE  CHANNEL 
in  Section  A1^2  by  the  following  set  of  relations: 
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Hi(8)  -  Hd(8)  H^,(8) 

H2(8)  -  Hh(8)  H^(8) 

H3(8)  -  Hu(8)  H^Cs) 

H^(8)  -  H^(8) 


In  addition,  noise  n(t)  ha8  been  added  in  the  tranamlaaion 
channel.  Thie  noise  is  assumed  to  have  a  spectral  density  of  ^  C^) 
which  may  be  white  compared'  to  the  system  bandwidth. 

Vnien  considering  the  signal  component,  X(t),  only,  the  double 
sampling  operation  can  be  replaced  by  a  single  sampling  operation  in  order 
to  achieve  even  greater  simplicity.  This  simplification,  which  can  be 
justified  rigorously,  can  also  be  obtained  heuristically  by  noting  that 
the  outputs  of  both  multipliers  are  sequences  of  impulses.  Hence,  the 
second  multiplier  and  the  transfer  function  between  it  and  the  first  can 
be  replaced  by  pure  attenuation  and  delay  only.  The  resulting  block  dia- 


FIG.  Al-3  SIMPLIFIED  BLOCK  DIAGRAM  FOR  SIOJAL  ONLY 

The  attenuation  factor  K  is  simply  the  magnitude  of  the  unit  im¬ 
pulse  response  of  H2(s)  H3  (s)  evaluated  at  the  receiver  sampling  times. 
Since  these  sampling  times  are  delayed  by  't'  seconds,  it  follows  that 

j«o 


since  this  is  just  the  Inverse  Laplace  transform  of  the  transfer  function 
evaluated  at  the  time 

It  will  be  shown  later  that  this  same  simplification  can  be  ex¬ 
tended  to  consider  crosstalk  by  using  different  values  for  ^  .  However, 

the  effects  of  noise  cavmot  be  considered  in  Figure  Al-3  and  it  is  necessary 
to  use  the  form  shown  in  Fi4,ure  Al-2  for  this  purpose. 
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An  assumption  vhlch  is  Implicit  in  the  use  of  Figure  Al-3  is 
that  the  Impulse  response  of  H2(s)H2(s)  Is  essentially  zero  at  a  time 
equal  to  the  sampling  Interval  T]^.  Since  small  crosstalk  requires  that 
this  Impulse  response  be  small  at  a  time  tj^  £l  ,  It  seems  clear  that 

N 

assumption  Is  extremely  well  Justified. 


Al-4  ERROR  ANALYSIS 


Errors  of  three  types  are  considered  in  this  analysis.  These 
are  designated  as  signal  error,  crosstalk,  and  noise.  The  signal  error 
arises  from  three  different  sources.  First,  there  is  an  error  due  to  the 
finite  bandwidth  of  the  data  filter  and  sampling  window.  Secondly,  there 
Is  the  aliasing  error  resulting  from  the  spectrum  folding  of  the  sampling 
process.  Finally,  there  is  a  reconstruction  error  resulting  from  imperfect 
Interpolation.  Since  all  of  these  signal  error  components  depend  upon  the 
channel  transfer  functions  and  the  sampling  Interval,  It  Is  neither  neces¬ 
sary  nor  desirable  to  treat  them  separately. 


For  purposes  of  error  analysis  It  will  be  assumed  that  the  sig¬ 
nal  X(t)  from  the  data  source  Is  a  random  function  with  a  spectral  density 
(S)  .  The  crosstalk  and  noise  are  also  random  functions  and  are 

assumed  to  be  statistically  Independent  of  each  other  and  of  the  signal. 


If  y(t)  is  the  output  from  any  channel  due  to  the  signal  input 
to  that  channel  only,  then  the  signal  error  may  be  defined  as 

s(t)  =  x(t  -  tp)  -  y(t)  (2) 

The  delay  time  to  is  introduced  in  order  that  delay  in  the  chan¬ 
nel  not  be  interpreted  as  an  error.  The  spectral  density  of  this  error 
may  be  written  in  terms  of  the  spectral  densities  of  x([t)  and  y(t)  and 
the  cross-spectral  densities.  Thus, 


$  CS) 


-  i 


i  (s')- 
^yy 


St. 

C 


-s-t« 


ts)-e 


</ 


^  (S) 


(3) 


The  spectral  densities  in  (3)  may  be  determined  from  the  sim¬ 
plified  block  diagram  of  Figure  Al-3.  The  only  item  of  special  concern 
is  to  correctly  account  for  the  effect  of  sampling.  This  may  be  done  with 
the  aid  of  Figure  Al-4  which  shows  the  multiplier  only. 
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FIG.  Al-4  MULTIPLIER  WHICH  IS  EQUIVALENT  TO  SAMPLING 


Since 

X2(t)  =  Xj(t)f^(t)  (4) 

the  spectral  density  of  X2(t)  is  obtained  by  convolution.  Thus, 


where 


d'JO 


(5) 


-  j.* 


^22^  =  Spectral  density  of  X2(t) 
^  (s)  =  spectral  density  of  Xj^(t) 
§  =  spectral  density  of  fiCt) 

TT  ^ 


Since  the  sampling  function  fj^Ct)  is  a  periodic  sequence  of  unit 
impulses,  its  spectral  density  can  be  shown  to  be 


(6) 


nr-'O© 


Substituting  this  function  into  (5)  yields 


-hi 

'  n-  -*0  ' 

Furthermore,  it  can  be  shown  that  cross-spectral  density  is 


* 

21 


(8) 


where  the  asteisk  denotes  the  complex  conjugate. 
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Using  (7),  (8)  and  Figure  Al-3,  the  signal  error  spectral  den 
sity  can  be  written  as 

-  —  e  H,<a^  H^(s)  ^  e  H,C'S)  H^(-s)  i 

The  error  due  to  crosstalk  can  be  obtained  by  extending  Figure 
Al-3  to  include  all  channels.  This  extension  is  shown  in  Figure 4l-S 
The  major  difference 


1 
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iH.Cs'i 

— {§)— 
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« 

^,Ct+(a-0t,) 

H^ii) 
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-sr 


FIG.  Al-5  BLOCK  DIAGRAM  FOR  CROSSTALK  ANALYSIS 


among  the  various  channels  is  the  attenuation  factor,  which,  for  the  kth 
channel ,  is 

,  f  CAI  (*) 

Hj  (s)  e  ds  (10) 

If  y(t),  in  Figure  Al-5,  is  considered  to  be  the  output  due  to 
X]^(t)  and  c(t)  is  the  crosstalk  output  due  to  all  other  Inputs,  then  the 
spectral  density  of  the  crosstalk  can  be  written  as 
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irz  »  =  .« 

\^ere  spectral  density  of  Xjj(t). 

The  error  due  to  noise  can  be  obtained  from  the  block  diagram 
of  Figure  Al-6. 

4(-fc) 


Tfltt) 


FIG.  Al-6  BLOCK  DIAGRAM  FOR  NOISE  ANALYSIS 

The  similarity  of  this  block  diagram  to  Figure  Al-3  makes  it  possible 
to  write  the  spectral  density  of  no(t)  by  analogy.  Thus, 


rfi.c.- 


,  -in  A 

I'T 


(12) 


Usually  the  noise  n(t)  and  the  transfer  function  H3(s)  will 
both  have  a  bandwidth  which  is  large  compared  to  the  sampling  frequency. 
In  this  case,  the  summation  in  (12)  can  be  closely  approximated  by  an 
integral.  If  this  is  done,  the  output  noise  spectral  density  can  be 
written  as 

Jeo 


ent, 

mean 


ds 


(13) 


Since  the  various  sources  of  error  are  statistically  independ- 
their  spectral  densities  can  be  added  directly.  Thus,  the  total 
squ.are  error  can  be  obtained  from 
joC 


I 


(14) 
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Upon  combining  some  terms  and  simplifying  the  notation,  this 
mean  square  error  can  be  written  as 


j«0 


-i-e 


-riff  k' 


where 


%]) 


c^S 


(• 


-4  — 

^  >"3  4.0 


It  may  be  noted  that  is  the  mean  square  value  of  noise  into 
the  receiver  sampler. 

Another  quantity  which  will  be  needed  in  the  remaining  analysis 
is  the  mean  square  value  of  the  signal  into  the  receiver  sampler.  This 
mean  square  value  is  given  by 

j*® 


s-  ^ 

3  2fT3 


'"3^1 


'SI  cis 


(16) 


where 


yaO 

Q  B  -  J  Hj  Ci)  Wj(rO] 


is  the  gain-bandwidth  product  of 


Al-5  ERROR  MINIMIZATION 

It  is  apparent  from  (15)  that  there  are  many  different  ways  in 
which  the  error  could  be  minimized.  For  purposes  of  the  present  analysis 
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'fv 


.1^' 


it  will  be  assumed  that  Hjt'Si,  HjC*!  ,  and  TJ  are  fixed  and  that  HiC®) 
and  can  be  adjusted  for  minimum  error.  It  will  also  be  assumed 
that  all  channels  are  identical  so  that 

N,  is)  (s)  =  HjCs') 

The  error  minimization  must  also  be  carried  out  with  a  constraint 
on  the  average  signal  power.  Without  this  constraint  the  effects  of  noise 
can  be  made  negligible  by  simply  increasing  the  signal  power.  Thus,  it  is 
necessary  to  minimize  the  function 


I  -  e'  4  XS. 


(17) 


where  X  is  the  Lagrange  multiplier  vdiich  can  be  adjusted  to  satisfy  the 
signal  power  constraint. 

The  notation  can  be  further  simplified  by  suppressing  the  inte¬ 
gration  variable  S  and  letting 


H.M  = 

n* 


etc. 


Furthermore,  let 


k:  - 1 


) 


The  integral  I  may  now  be  written  as 

\  n 

4  133.  U.uf  ]  ds  C18) 

Ti 


Further  work  is  simplified  if  the  optimum  value  of  tg  is  chosen. 
It  can  be  shown  straightforwardly  that 


IL 

St, 


-  o 


when  tg  =  T  .  This  value  of  tg  will  be  assumed  throughout  the  rest 
of  the  analysis. 

The  next  step  is  to  take  a  first  variation  of  I  with  respect 
to  H^.  Setting  5 1=  0  and  eliminating  the  redundant  conjugate  terms 
yields 

-  'in  "  J 


from  which 


(20) 


Before  taking  a  variation  with  respect  to  it  is  convenient 
to  mal^e  a  change  in  variable  in  the  portion  of  the  integral  containing 
Htn  •  The  resulting  integral  contains  a  summation  over  |-j^^  , 

Hence , 


1- 


371 


it  to  0, 

J-O 


Ti 

After  taking  the  first  variation  with  respect  to  and  equating 

-r  n*  .  ><56 


from  which 


H,1  %  I K,,  ^  ^  $.1  .  ^  iX 


(22) 


Equations  (21)  and  (22)  specify  the  conditions  that  must  be 
satisfied  in  order  for  the  coca!  mean  square  error  to  be  a  minimum. 
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Al-6  IDEAL  N0N-REALI2ABLE  TRANSFER  FUNCTIONS 


The  simultaneous  equations  which  must  be  satisfied  by  and 
H4  may  be  written  as 


and 


k[<1  I 


(23) 


(24) 


A  general  solution  to  this  pair  of  equations  would  appear  to 
be  quite  formidable.  If  it  is  assumed,  however,  that  and  H4  are  band- 
limited  so  that 


Hi 

II 

n  =  0 

=  0 

n  y  0 

H4  Hin 

=  H1H4 

n  =  0 

=  0 

n  4  0 

then  a  solution  is  readily  obtained.  This  assumption  may  be  justified, 
a  posteriori .  by  noting  that  the  resulting  solutions  are,  in  fact,  band- 
limited. 


The  solutions  obtained  by  this  procedure  are 


Al-13 


which  exist  when  -V  and  are  zero  otherwise.  The  value  of 

^  is  adjusted  by  putting  (25)  into  (16)  with  the  desired  value  for 
S3.  This  may  then  be  solved  for  ^  . 

The  solutions  given  by  (25)  and  (26)  are  in  general  not  physi¬ 
cally  realizable  because  the  resulting  transfer  functions  may  have  RHP 
poles  or  complex  poles  not  occurring  in  conjugate  pairs.  In  addition, 
the  resulting  transfer  functions  may  not  be  rational.  Nevertheless, 
these  transfer  functions  are  useful  in  specifying  the  irreducible  error 
in  the  system.  This  irreducible  error  is  obtained  by  substituting  (25) 
and  (26)  into  the  expression  for  the  mean  square  error  as  given  by  (15) 
The  general  result  is  too  complicated  to  be  useful  but  specific  cases 
can  be  worked  out. 


Al-7  OPTIMUM  REALIZABLE  TRANSFER  FUNCTIONS 

The  optimum  realizable  transfer  functions  are  those  which  pro¬ 
duce  minimum  error  and  have  all  poles  in  the  LHP.  In  general,  these 
transfer  functions  would  require  an  infinite  number  of  poles  in  order  to 
achieve  the  bandlimited  condition  and  can  only  be  approximated  by  practi¬ 
cal  networks.  This  approximation  will  not  be  discussed  in  this  report. 

In  accordance  with  the  usual  procedures  for  obtaining  realizable 
functions,  it  is  necessary  to  introduce  a  shaping  filter  which  yields  a 
constant  spectral  density  at  its  output.  This  shaping  filter  transfer 
function  is  then  extracted  from  the  ideal  transfer  function  and  the  re¬ 
mainder  is  made  realizable  by  retaining  only  its  LHP  poles. 

Considering  the  data  filter  first,  the  required  shaping  filter 
has  a  transfer  function  Hg(s)  defined  by 


or 


HiO'jUsC-i')  - 


(27) 


I  ^  T  *  ^ 

Hi  Hi  = 

Thu  remainder  ol  ,  after  extracting  Hg ,  is  designated  as 
H]Ls  is  given  by 


(28) 
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from 


The  physically  realizable  part  of  is  and  is  obtained 


+  Im-I 


(29) 


by  making  a  partial  fraction  expansion  of  if  it  is  rational,  or  by 
performing  the  equivalent  time-domain  separation  if  it  is  not.  In 
either  case  [His1+  contains  only  LHP  poles  and  [H]^g^  .  contains  only  RHP 
poles.  The  resulting  overall  realizable  transfer  function,  H^r  is  now 
given  by 


(30) 


Upon  using  (25)  and  defining  some  new  terms,  the  realizable  fil¬ 
ter  becomes 


where 


K]± 

4>. 


(31) 


Aa  e> 


It  should  be  noted  that  Hj^j^  contains  the  transfer  function  which 
corresponds  to  the  sampling  window  as  well  as  the  data  filter.  This  im¬ 
plies  that  the  data  filter  Hjj  is  given  by 


- 


U  W 


(32) 


The  only  restriction  on  the  sampling  window,  therefore,  is  that 
Hy  must  not  have  transmission  zeros  within  the  passband  of  Hj^j^.  This 
condition  is  almost  always  satisfied. 

Since  the  input  to  H^  is  a  sequence  of  impulses,  no  shaping 
filter  is  needed.  Hence  the  realizable  filter  transfer  function  is 


"  W+  ■ 


(33) 
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where 


Al-8  OPTIMUM  SAMPLING  RATE 


It  was  shown  in  the  previous  sections  that  the  optimum  data 
filter  and  interpolation  filter  were  bandlimited.  Hence,  any  sampling 
rate  greater  than  that  determined  by  the  filter  bandwidth  will  not  be 
effective  in  reducing  signal  error.  However,  increased  sampling  rates 
will  increase  the  error  due  to  noise  because  the  video  bandwidth  must 
increase  (for  a  fixed  number  of  channels). 

It  may  be  concluded  from  the  above  that  the  optimum  sampling 
rate  is  that  determined  by  the  filter  bandwidth.  The  filter  bandwidth 
is  that  which  makes  in  (31)  and  (32)  identically  zero.  Thus,  if 

- 

then 


Al-9  OPTIMUM  SAMPLING  DELAY  AND  VIDEO  FILTERS 


>  68  N: 


(34) 


A  set  of  attenuation  factors  was  defined  by  (10)  as 
iU)  (t^  S(T+(4-Oi,') 

H,  e  ds 


2T1 


r 

_  Hi  CM 

i  --jo* 


(35) 


These  factors  may  be  interpreted  in  the  time  domain  by  sketch¬ 
ing  the  impulse  response  corresponding  to  H2(s)  H3(s)  under  the  assump¬ 
tion  that  the.se  are  the  same  for  all  channels.  The  Kj^  are  then  simply 
samples  of  this  impulse  response,  as  shown  in  Figure  Al-7.  It  is  clear 
from  the  mean  square 
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error  as  given  by  (15),  or  simply  from  the  physical  interpretation, 
that  Ki  should  be  as  large  as  possible  and  all  of  the  others  should  be 
as  small  as  possible. 

The  optimum  value  of  sampling  delay,  ,  should  correspond 
to  the  location  of  the  maximum  of  the  impulse  response.  The  following 
zeros  of  the  Impulse  response  should  coincide  with  ^  ~C  i-2'ti  j  -  .  . 

etc.  Hence,  if  the  video  filters  have  been  selected  so. as  to  minimize 
crosstalk,  the  optimum  value  of  'f  is  that  which  satisfies  the  equation 


d 

d  t 


-  O 


(36) 
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APPENDIX  2 


POWER  SPECTRUM  OF  SAMPLER  OUTPUT 


The  sampling  operation  is  considered  as  equivalent  to  taking 
the  product  of  fd(t)  +  N(tX] •  s(t)  =  g(t) ,  where  s(t)  is  the  sampling 
function  shown  in  Figure  3-2.  The  Fourier  transform  of  8(t)  is  deter¬ 
mined  below: 
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The  power  spectrum  of  s(t)  is  given  by 
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where  the  bar  denotes  an  ensemble  average. 
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Note  that  if  A  =  tj^  and  tj^ — >-0,  then  PgCf)  =  Tg^  ^  0(f  -  Tg  j 

which  is  the  same  result  obtained  if  the  sampling  function  is  a  train 
of  delta  functions  in  the  time  domain. 
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The  power  spectrum  of  the  sampler  output  signal  is  determined 
by  convolvin'’  the  power  spectrum  of  the  sampling  function  with  the  power 
spectrum  of  the  sampler  input.  Let  h(t)  =  d(t)  +N(t) 
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LIST  OF  REPORTS  AND  CONFERENCES  ON  DATA  SPECTRA  SURVEY 


1)  "Basic  Data  Compilation  -  Short  Range  Flight  Test  Operation 
No.  5",  22  February  1961,  Aeronutronic  Document  18640-N 
(Confidential) 

2)  Aliasing  Study  -  Review  of  Report 

Title:  Study  and  Evaluation  of  In-Flight  Guided  Missile  Environ¬ 
mental  Data 

Date:  (1)  Interim  Report  #1,  Report  #F161g-l  (63248),  13  June  1958 
(2)  Final  Report,  Report  #F161-72809,  ASTIA  AD  306  231, 

15  December  1958 

Reports  Confidential  but  title  unclassified 
Source:  Fairchild  Engine  &  Airplane  Corporation 
Applied  Science  Division 
1334  North  Henry  Street 
Alexandria,  Virginia 
Prepared  for:  Navy  Bureau  of  Aeronautics 
Contract  No.  NOar  58-539C 
This  report  includes  an  extensive  bibliography. 

3)  Wiancko  Engr.  Rpt.  #WC2329-F  dated  May  1959  "S'.condary  Isolation 
System  for  the  High  Speed  Track  Edwards  Air  Force  Base" 

4)  Northrop  Aircraft  Inc.  Report  NAI-57-586  "Summary  of  Launch  and  Flight 
Acoustic  and  Vibration  Data  N-69D  Missiles",  dated  12  March  1959 

5)  "Basic  Data  Compilation  Fire  Control  Static  Test  I-S-2", 

19  April  1961,  Aeronutronic  D.C.  20025-N  (Secret) 

6)  Spectrum  Survey  Trip  Report  -  JPL  -  Vibration  Analysis  Lab. 
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BIBLIOGRAPHIES  OBTAINED  FOR  DATA  SPECTRA  SURVEY 


Bibliography  of  Electronics  Reports  for  M.I.T. 

Research  Laboratory 
February  1 ,  1958 

Approximately  338  reports  referenced 

A  Bibliography  of  Telemetry 

IRE  Transactions  on  Telemetry  and  Remote  Control 
June ,  1958 

Approximately  533  books,  reports  and  papers  referenced 
A  Report  Bibliography  on  Telemetering  and  Telemeter  Systems 
AD  55  316 

Bibliography  Secret 
Includes  280  references 
Compiled  by  H.  L.  Conn,  May  1954 

A  Report  Bibliography  on  Telemetering  and  Telemeter  Systems 

AD  30  076 

Bibliography  Secret 

Includes  approximately  250  references 
Compiled  by  H.  L.  Conn,  May  1954 
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ASTIA  BIBLIOGRAPHIES 
CODE:  6065-812-I/S 


1.  A  Report  Bibliography  on  Phase  Detectors 

ARB-6440 
7  secret  cards 
14  confidential  cards 
54  unclassified  cards 

2.  Bibliography  on  "Errors"  Communication  Systems 

ARB-6447 
39  secret  cards 
34  confidential  cards 
75  unclassified  cards 

3.  A  Report  Bibliography  on  Data  Processing  Systems 

ARB-6445 
219  secret  cards 

112  confidential  cards 
604  unclassified  cards 

4.  Report  Bibliography  on  Data  Transmission  Systems 

ARB-6446 

113  secret  cards 

234  confidential  cards 
455  unclassified  cards 

5.  A  Report  Bibliography  on  Phase  Measurement 

ARB-6443 
2  secret  cards 
6  confidential  cards 
Approx.  65  unclassified  cards 

6.  A  Report  Bibliography  on  Phase  Modulation 

ARB-6444 
13  secret  cards 
6  confidential  cards 

36  upfclasslfied  cards 

7.  Report  Bibliography  on  "Low  Pass  Filters" 

ARB-6441 

6  secret  cards 

7  confidential  cards 
Approx.  60  unclassified  cards 
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8.  A  Report  Bibliography  on  Telemetering  Data 

ARB~eiA2 
16  secret  cards 
21  confidential  cards 
90  unclassified  cards 

9.  Atlantic  Missile  Range 

U.  S.  Air  Force  Missile  Test  Center 
Test  Activity  Report 
September,  1960 

(Title  unclassified  but  Report  Secret) 

AD  319  619 

10.  A  Report  Bibliography  on  Radio  and  Communication  System 
Interference 

ARB-6681 
ARB-6682 
ARB-6680 
105  Secret  cards 
70  confidential  cards 
520  unclassified  cards 
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